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Abstract:  After  much  previous  research,  many  controversial  issues  re¬ 
lated  to  well  developed  design  recommendations  and  well  accepted  ana¬ 
lytical  procedures  for  infilled  frames  are  still  unresolved.  The  main  diffi¬ 
culty  in  evaluating  the  performance  of  bearing  and  shear  wall  structures  is 
determining  the  type  of  interaction  between  the  infill  and  the  frame,  which 
has  a  major  impact  on  the  structural  behavior  and  load-resisting  mecha¬ 
nism.  A  simplified  model  that  captures  the  characteristics  of  infill  masonry 
would  be  valuable,  but  developing  one  requires  understanding  of  ma¬ 
sonry-infilled  frame  behavior  in  much  more  detail  than  strut  or  beam  be¬ 
havior.  Numerical  simulations  using  well  calibrated  constitutive  models 
are  needed  for  parametric  studies  to  facilitate  the  development  and  cali¬ 
bration  of  a  simple,  accurate  infill  masonry  model.  The  objective  of  this 
study  was  to  identity  suitable  numerical  constitutive  models  and  demon¬ 
strate  their  capabilities.  Two  were  studied  here:  a  cohesive  interface  model 
to  simulate  the  behavior  of  mortar  joints  between  masonry  units  as  well  as 
the  behavior  of  the  frame-to-panel  interface;  and  a  smeared  crack  finite 
element  formulation.  The  interface  model  was  able  to  account  for  the 
shearing,  residual  shear  strength,  and  opening  and  closing  of  joints  under 
cyclic  shear  loads  in  simple  combinations  of  concrete  blocks  and  mortar 
joints. 
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Notations 

Cs  =  parameter  controlling  the  shear  stress  contribution  to  failure,  DIANA  interface 
model 

E  =  modulus  of  elasticity,  psi,  DIANA  smeared  crack  model 
Gfc  =  fracture  energy  in  compression,  psi-in.,  DIANA  models 
Gf1  =  first  mode  fracture  energy,  psi-in. 

Gf11  =  shear  mode  fracture  energy,  psi-in. 

Knn  =  normal  numerical  elastic  stiffness  parameter,  psi/in. 

Kss  =  tangential  numerical  elastic  stiffness  parameter,  psi/in,  DIANA  interface  model 
c0  =  initial  cohesion,  psi,  DIANA  interface  model 
f  c  =  compressive  strength,  psi,  DIANA  models 
f  cm  =  compressive  strength  of  masonry  mortar  cubes,  psi 

fcu  =  compressive  strength  of  masonry  units,  psi  (with  respect  to  net  cross-sectional  area) 

f  m  =  compressive  strength  of  masonry  prisms,  psi  (with  respect  to  net  cross-sectional 
area) 

ft  =  joint  tensile  strength,  psi,  DIANA  models 

\p  =  dilatancy  angle,  radian,  DIANA  interface  model 

P  =  shear  retention  factor,  DIANA  smeared  crack  model 

5  =  dilatancy  shear  slip  degradation  coefficient,  DIANA  interface  model 

eu  =  strain  at  maximum  strength  for  axial  testing  of  masonry  prisms  (in./in.) 

4»  =  internal  friction  angle,  radian,  DIANA  interface  model 

4>i  =  initial  internal  friction  angle,  radian,  DIANA  interface  model 

4>r  =  residual  internal  friction  angle,  radian,  DIANA  interface  model 

kp  =  norm  of  plastic  strain  associated  with  peak  compressive  strength,  DIANA  interface 
model 

v  =  Poisson’s  ratio,  DIANA  smeared  crack  model 


On  =  confining  stress  over  which  the  dilatancy  will  be  zero,  psi,  DIANA  interface  model 
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1  Introduction 

Background 

Masonry  walls  have  long  been  used  as  both  load-bearing  structural  ele¬ 
ments  and  architectural  non-structural  elements  in  single-  and  multi-story 
buildings.  Both  reinforced  and  unreinforced  masonry  have  been  used  in 
such  structures  on  military  installations.  In  some  cases,  masonry  walls  fill 
the  space  within  structural  frames  to  separate  interior  spaces.  In  other 
cases  load-bearing  masonry  walls  are  not  bounded  by  structural  frames, 
but  are  joined  with  the  structural  frame  using  ties  and  columns  to  provide 
better  structural  integrity  to  the  masonry. 

Masonry  infills  are  usually  treated  as  non-structural  elements  in  buildings, 
and  their  interaction  with  the  bounding  frame  is  often  ignored  in  design. 
Nevertheless,  infills  contribute  strength  to  a  structure  and  will  interact 
with  the  bounding  frame  when  the  structure  is  subjected  to  strong  lateral 
seismic  loads.  This  interaction  may  or  may  not  be  beneficial  to  the  per¬ 
formance  of  the  structure,  however,  and  it  has  been  a  topic  of  much  debate 
in  the  last  few  decades. 

Infill  walls  have  been  identified  as  a  contributing  factor  to  catastrophic 
structural  failures  in  earthquakes.  Frame/partial  infill  interaction  can 
cause  brittle  shear  failures  of  reinforced  concrete  columns  and  short- 
column  effect.  Furthermore,  infills  can  over-strengthen  the  upper  stories 
of  a  structure  and  result  in  a  “soft”  first  story,  which  is  highly  undesirable 
from  the  perspective  of  seismic  safety. 

In  spite  of  the  shortcomings  that  have  been  observed  in  this  type  of  con¬ 
struction,  there  is  strong  laboratory  and  field  evidence  that  masonry  infills 
can  improve  the  earthquake  resistance  of  a  frame  structure  if  they  are 
properly  designed. 

Literature  review 

The  behavior  of  masonry-infilled  steel  and  reinforced  concrete  frames  un¬ 
der  in-plane  lateral  loading  has  been  the  subject  of  many  investigations. 
Fiorato  et  al.  (1970)  tested  1/8-scale  non-ductile  reinforced  concrete 
frames  infilled  with  brick  masonry  under  monotonically  increasing  as  well 
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as  cyclic  lateral  loading.  Those  tests  were  followed  by  the  studies  of  Kling- 
ner  and  Bertero  (1976),  Bertero  and  Brokken  (1983),  Zarnic  and 
Tomazevic  (1985),  and  Schmidt  (1989).  More  recently,  single-story  rein¬ 
forced  concrete  frames  with  masonry  infills  were  studied  by  Mehrabi  et  al. 
(1994, 1997),  Angel  et  al.  (1994),  and  Al-Chaar  et  al.  (1998,  2002).  The  lat¬ 
ter  study  also  examined  the  behavior  of  masonry  infill  under  out-of-plane 
loading. 

Masonry  infilled  steel  frames  were  tested  by  Dkanasekar  et  al.  (1985), 
Dawe  and  Seah  (1989),  Mander  et  al.  (1993),  Mosalam  and  Paulino  (1997), 
and  Buonopane  and  White  (1997).  Flanagan  and  Bennett  (1999a)  investi¬ 
gated  steel  frames  with  structural  clay  tile  infills. 

All  studies  cited  above  have  shown  that  the  behavior  of  an  infilled  frame  is 
heavily  influenced  by  the  interaction  of  the  infill  with  its  bounding  frame. 
In  most  instances,  the  lateral  resistance  of  an  infilled  frame  is  not  equal  to 
a  simple  sum  of  the  resistance  if  its  components  because  frame/infill  in¬ 
teraction  can  alter  the  load-resisting  mechanisms  of  the  individual  com¬ 
ponents.  At  low  lateral  loading,  an  infilled  frame  acts  as  a  monolithic  load- 
resisting  system.  However,  as  loading  increases,  the  infill  tends  to  partially 
separate  from  the  bounding  frame  and  form  a  compression  stmt  mecha¬ 
nism,  as  observed  in  many  earlier  studies  (e.g.,  Stafford  Smith  1962). 
However,  the  compression  stmt  may  or  may  not  evolve  into  a  primary 
load-resistance  mechanism  for  the  structure,  depending  on  the  strength 
and  stiffness  properties  of  the  infill  with  respect  to  those  of  the  bounding 
frame. 

Dhanasekar  and  Page  (1986)  and  Liauw  and  Lo  (1988)  have  used  linear 
and  nonlinear  beam  elements  to  model  the  behavior  of  steel  frames,  and 
interface  elements  to  model  the  interaction  between  the  infill  and  the 
frame.  Dhanasekar  and  Page  used  a  nonlinear  orthotropic  model  to  simu¬ 
late  the  behavior  of  brick  infills,  and  Liauw  and  Lo  used  a  simple  smeared 
crack  model  to  simulate  the  behavior  of  micro-concrete  infills.  Schmidt 
(1989)  used  smeared  crack  elements  to  model  both  reinforced  concrete 
frames  and  brick  infills.  In  all  those  analyses,  infill  panels  have  been  mod¬ 
eled  as  a  homogenous  material  before  fracture,  and  the  effects  of  mortar 
joints  have  been  smeared  out.  While  smeared  crack  models  have  been 
used  to  analyze  reinforced  concrete  structures  where  cracks  are  often  dif¬ 
fuse,  they  suffer  from  a  major  deficiency  in  that  they  cannot  properly  cap¬ 
ture  diagonal  shear  cracks  and  the  shear  sliding  of  cracked  concrete  or  ma- 
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sonry  mortar  joints.  Those  deficiencies  are  inherent  in  the  kinematic  con¬ 
straints  related  to  trying  to  account  for  cracking  in  a  continuum  (i.e.,  a 
homogeneous  material).  The  use  of  smeared  crack  elements  alone  will  lead 
to  non-conservative  design  results  for  unreinforced  masonry  infill  because 
the  masonry  has  natural  planes  of  weakness  introduced  by  mortar  joints. 
The  most  realistic  approach  to  accounting  for  those  weak  planes  is  to  in¬ 
corporate  interface  elements  into  the  model. 

A  number  of  plasticity-based  continuous-interface  models  have  been  de¬ 
veloped  to  model  the  tension  and  shear  behavior  of  masonry  mortar  joints 
(Rots  1991;  Lotfi  and  Shing  1994;  Lourenco  1996).  Those  models  account 
for  the  interaction  between  normal  compression  and  shear  as  well  as  the 
shear  dilatation  often  observed  in  experiments.  Mehrabi  and  Shing  (1997) 
have  developed  an  interface  model  for  analyzing  masonry  infills  that  ac¬ 
counts  for  the  increase  of  contact  stress  due  to  joint  closing,  the  geometric 
shear  dilatation,  and  the  plastic  compaction  of  a  mortar  joint.  The  failure 
surface  of  the  model  is  based  on  a  hyperbolic  function  proposed  by  Lotfi 
and  Shing  (1994),  and  is  capable  of  modeling  damage  accumulation  at 
mortar  joints  under  increasing  displacement  and  cyclic  loading.  This  is  re¬ 
flected  by  shear  strength  reduction  and  mortar  compaction  (loss  of  mate¬ 
rial)  at  interfaces.  The  model  has  been  used  to  analyze  the  infilled  frames 
tested  by  Mehrabi  et  al.  (1994). 

Objective 

The  objective  of  this  study  was  to  carry  out  numerical  simulation  and  pa¬ 
rametric  studies  and  compare  them  with  the  results  of  a  previous  experi¬ 
mental  investigation.  The  investigation  also  attempted  to  determine  the 
type  of  interaction  between  the  infill  and  the  frame,  which  has  a  major  im¬ 
pact  on  the  structural  behavior  and  load-resisting  mechanism,  and  would 
capture  characteristics  of  reinforced  concrete  infilled  with  masonry  in 
much  more  detail  than  a  strut  or  beam. 

Approach 

A  review  of  the  literature  was  performed  to  summarize  the  state  of  the 
technical  domain  in  terms  of  published  experimental  and  analytical  re¬ 
search  results. 

Advanced  finite  element  modeling  approaches  were  selected  for  detailed 
analysis  of  the  nonlinear  behavior  of  infilled  frames.  The  capabilities  of 
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various  numerical  models  are  described.  Also,  to  facilitate  the  use  of  these 
models  by  structures  researchers  and  designers,  a  commercial  finite  ele¬ 
ment  program  with  similar  analytical  capabilities  was  identified  and 
tested. 

The  finite  element  models  used  in  this  study  were  (l)  a  cohesive  interface 
model  to  simulate  the  behavior  of  mortar  joints  between  masonry  units 
and  the  behavior  of  the  frame-to-panel  interface,  and  (2)  a  smeared  crack 
finite  element  formulation  to  model  the  concrete  in  R/C  frames  and  ma¬ 
sonry  units.  The  interface  model  was  used  to  analyze  simple  combinations 
of  concrete  blocks  and  mortar  joints  because  it  can  account  for  the  shear¬ 
ing,  residual  shear  strength,  and  opening  and  closing  of  joints  under  cyclic 
shear  loads.  The  analytical  models  are  compared  with  experimental  mod¬ 
els  used  in  previous  investigations  to  verily  results  and  validate  used  ana¬ 
lytical  parameters.  The  numerical  models  were  compared  with  previous 
experimental  results  to  verify  the  adequacy  of  assumptions  employed  in 
the  numerical  models  and  to  compare  results. 

Mode  of  technology  transfer 

The  constitutive  models  employed  in  this  investigation  can  be  applied  by 
the  qualified  structural  engineer  to  the  analysis  of  critical  and  complex 
masonry-infilled  R/C  frames  using  the  DIANA  finite  element  modeling 
and  simulation  program.  The  results  may  be  expected  to  have  an  accept¬ 
able  level  of  accuracy,  but  must  be  considered  within  the  context  of  the 
several  caveats  discussed  in  the  text. 
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2  Selection  of  Analytical  Methods  and  Tools 

Failure  mechanisms  of  infilled  panels  under  in-plane  lateral  loading 

Studies  have  shown  that  infilled  frames  can  develop  a  number  of  possible 
failure  mechanisms  depending  on  the  strength  and  stiffness  of  the  bound¬ 
ing  frames  with  respect  to  the  strength  and  stiffness  of  the  infills  and  the 
geometric  configuration  of  the  framing  system.  Most  analytical  models 
proposed  today  focus  on  one  type  of  mechanism  or  the  other,  but  neither 
is  universally  applicable  to  all  infilled  structures.  However,  recent  research 
has  shed  more  light  on  the  behavior  of  infilled  frames  and  has  resulted  in 
the  development  of  advanced  analytic  tools.  Out-of-plane  behavior  of  in¬ 
filled  frames  has  been  studied  in  addition  to  the  in-plane  response.  Classi¬ 
cal  diagonal  strut  models  have  been  subjected  to  more  thorough  evalua¬ 
tions  with  new  experimental  data,  and  various  limit  analysis  methods  have 
been  developed  to  account  for  the  different  load-resisting  mechanisms  of 
infilled  frames. 

Based  on  experimental  observations,  one  can  identify  five  main  failure 
mechanisms  of  infilled  frames.  They  are  illustrated  in  Figure  l  (Mehrabi  et 
al.  1994).  While  most  of  the  studies  to  date  have  focused  on  unreinforced 
masonry  panels,  Klingner  and  Bertero  (1976)  and  Bertero  and  Brokken 
(1983)  investigated  the  behavior  of  engineered  infilled  frames.  Both  inves¬ 
tigations  tested  one-third  scale,  three-story  reinforced  concrete  frames  in¬ 
filled  with  fully  grouted  concrete  masonry  that  had  both  horizontal  and 
vertical  reinforcement.  The  infill  panels  were  securely  tied  to  the  bounding 
frames.  Both  studies  demonstrated  that  properly  engineered  infilled 
frames  can  provide  superior  performance,  in  terms  of  strength,  stiffness, 
and  energy  dissipation,  as  compared  with  a  bare  frame.  However,  rein¬ 
forced  infills  are  not  common,  and  an  over-reinforced  infill  may  create  a 
risk  of  brittle  shear  failure  in  the  bounding  reinforced  concrete  columns. 
Studies  by  Mehrabi  et  al.  (1994, 1996)  have  shown  that  relatively  weak  un¬ 
reinforced  masonry  infills  can  significantly  enhance  the  stiffness  and 
strength  of  a  non-ductile  reinforced  concrete  frame  without  jeopardizing 
ductility. 
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In  summary,  the  failure  mechanism  and  load-resistance  of  an  infilled 
frame  depend  very  much  on  the  strength  and  stiffness  of  the  infill  with  re¬ 
spect  to  the  bounding  frame.  It  is  evident  that  the  strength  of  the  mortar 
joints  is  also  a  dominant  factor.  A  relatively  weak  infill  is  most  desirable. 
Studies  have  shown  that  infill  panels  can  significantly  enhance  the  per¬ 
formance  of  a  bare  frame  under  earthquake  loads  provided  the  short  col¬ 
umn  phenomenon  and  the  brittle  shear  behavior  of  the  columns  can  be 
prohibited  (Bertero  and  Brokken  1983;  Mehrabi  et  al.  1994). 

An  accurate,  computer-assisted  analytical  method  for  assessing  alternate 
infill  panel  designs  would  be  much  more  affordable  than  the  construction 
and  testing  of  individual  scaled  physical  models  for  a  specific  application. 
Owing  to  the  compound  complexities  of  the  modeling  and  simulation 
problem,  however,  two  fundamental  questions  first  needed  to  be  an¬ 
swered: 

•  What  is  the  most  suitable  analytical  approach  to  the  problem? 

•  What  available  numerical  models  best  fit  the  requirements  of  the  se¬ 
lected  analytical  approach  and  the  physical  characteristics  of  masonry- 
infilled  panels? 

Available  analytical  modeling  approaches 

The  candidate  analytical  methods  fall  into  three  categories: 

1.  simple  modeling 

2.  limit  analysis 
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3.  constitutive  micro-  or  macro-mechanical  models  for  finite  element  analy¬ 
sis. 

In  the  simple  modeling  method,  infill  panels  are  usually  simplified  as  one 
or  more  distinct  elements  in  an  attempt  to  capture  the  important  behav¬ 
ioral  parameters  of  an  actual  infilled  frame  structure.  Limit  analysis  meth¬ 
ods  are  based  on  a  lower-bound  solution  for  lateral  strength  of  an  infilled 
frame  based  on  a  series  of  assumed  failure  mechanisms.  Constitutive  finite 
element  modeling  methods  attempt  to  simulate  the  “actual  behavior”  of 
infill  panels  by  dividing  the  panels  into  finite  elements  and  implementing 
constitutive  models  for  these  elements  that  could  capture  the  mechanical 
characteristics  of  the  material  in  infill  panels.  In  this  paper,  only  constitu¬ 
tive  models  for  finite  element  analysis  are  reviewed. 

Evaluation  of  constitutive  model  alternatives 

Finite  element  analysis  may  be  the  most  powerful  tool  for  evaluating  the 
behavior  of  infilled  frames.  However,  the  reliability  of  the  finite  element 
method  depends  on  the  appropriateness  of  the  model  and  elements  used. 
A  properly  developed  finite  element  model  can  account  for  all  possible 
failure  mechanisms  of  infilled  frames,  but  improper  use  of  such  models 
can  produce  non-conservative  results.  Therefore,  extreme  care  must  be 
exercised  when  performing  a  finite  element  analysis. 

Many  different  approaches  have  been  used  for  finite  element  modeling  of 
masonry  infills.  Recent  approaches  have  recognized  the  inelastic  and  plas¬ 
tic  behavior  of  masonry  from  the  initial  stage  of  loading  and  the  effect  that 
damage  accumulation  has  on  the  modeling  of  the  masonry  infills.  The  re¬ 
view  presented  here  discusses  models  that  consider  nonlinear,  plastic  be¬ 
havior  and  damage  effects  resulting  from  masonry  infill  as  an  isotropic  or 
orthotropic  brittle  or  quasi-brittle  material  and/or  brittle  mortar  joints. 

Simpler  models  have  considered  the  masonry  infill  to  be  a  homogeneous 
isotropic  continuum  medium  with  no  tensile  strength  (Del  Piero  1989)  or 
by  applying  the  smeared  crack  approach.  However,  more  detailed  consti¬ 
tutive  models  must  be  considered  to  describe  the  unique  characteristics  of 
masonry  infills,  namely  a  highly  orthotropic  nature  due  to  the  presence  of 
mortar  joints  and  shear  strength  degradation  mostly  focused  along  mortar 
joints.  The  first  task  in  incorporating  more  detail  is  to  introduce  appropri¬ 
ate  constitutive  models  for  both  the  masonry  units  and  mortar  joints.  To 
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this  end,  three  different  approaches  can  be  considered  based  on  different 
levels  of  complexity  (Gambarotta  and  Lagomarsino  1997): 

1.  modeling  masonry  units  and  mortar  joints  as  continuum  elements 

2.  continuum  modeling  of  combination  units  and  joints  using  a  homogeniza¬ 
tion  procedure 

3.  introduction  of  interface  elements. 

A  more  effective  approach,  however,  considers  mortar  joints  —  the  weak¬ 
est  elements  of  unreinforced  masonry  —  as  interface  elements  (Lotfi  and 
Shing  1994;  Mehrabi  and  Shing  1997).  The  constitutive  modeling  of  inter¬ 
faces  has  improved  over  time  to  accommodate  the  orthotropic  behavior  of 
masonry  and  model  separation,  shear  degradation,  and  closing  due  to 
unloading.  The  interfaces  connect  masonry  units.  Loading  caused  by 
unloading  can  be  modeled  using  available  smeared  crack  models  with  or 
without  damage  consideration.  The  presence  of  interfaces  overcomes  the 
shear  locking  and  mesh  dependency  often  associated  with  the  use  of 
smeared-crack  models  alone. 

There  exist  various  approaches  to  the  constitutive  modeling  of  brittle  and 
quasi-brittle  materials.  Many  use  continuum  damage  modeling.  Some  ad¬ 
dress  the  modeling  of  concrete-type  material  and  not  necessarily  compos¬ 
ite  masonry  material.  However,  those  models  form  the  basis  for  the 
smeared  crack  models  that  can  be  effectively  used  for  modeling  of  ma¬ 
sonry  units  and  R/C  frames.  Several  other  models  directly  address  the 
modeling  of  masonry  in  general  and  some  masonry  infilled  frames  in  par¬ 
ticular.  However,  a  review  of  continuum  damage  modeling  is  beyond  the 
scope  of  this  report.  It  should  be  noted  that  for  the  case  of  unreinforced 
concrete  masonry,  infills  are  not  necessarily  as  sensitive  to  damage  in  ma¬ 
sonry  units  as  to  damage  in  mortar  joints.  Nevertheless,  the  majority  of 
smeared  crack  models  introduced  in  recent  years  consider  one  or  another 
type  of  damage  propagation  and  permanent  mechanical  property  varia¬ 
tion.  Some  of  the  interface  models  discussed  below  effectively  consider  the 
damage  initiation  and  progress  in  mortar  joints,  and  therefore  offer  the 
best  candidates  from  which  to  select  a  best-fit  model  for  unreinforced  con¬ 
crete  masonry  infilled  frames. 

Candidate  models  for  concrete,  masonry,  and  interface  elements 

This  section  reviews  in  more  detail  work  presented  by  several  researchers 
addressing  finite  element  modeling  of  brittle  and  quasi-brittle  materials 
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including  masonry  infills.  The  continuum  models  presented  here  mostly 
consider  the  unique  characteristics  of  such  material,  e.g.,  cracking  and 
damage,  smeared  over  a  finite  element. 

Mosalam  and  Paulino  (1997)  employed  fixed  smeared-crack  concepts  with 
strain  decomposition  to  carry  out  nonlinear  finite  element  simulations  of 
either  unreinforced  or  reinforced  concrete  structures.  To  alleviate  spurious 
mesh  dependence,  researchers  adapted  the  finite  element  mesh  to  the  pre¬ 
sent  level  of  damage  (cracking)  in  space  and  time.  Wiliam  and  Warnke 
(1974)  described  different  models  for  the  failure  surface  and  the  constitu¬ 
tive  behavior  of  concrete  under  triaxial  stress.  Dhanasekar  et  al.  (1985)  de¬ 
rived  simple  nonlinear  stress/strain  relations  for  brick  masonry  con¬ 
structed  with  solid  pressed  bricks  from  the  results  of  a  large  number  of 
biaxial  tests  on  square  panels  using  various  angles  of  the  bed  joint  in  rela¬ 
tion  to  the  principal  stress  axes.  Macroscopic  elastic  and  nonlinear 
stress/strain  relations  were  determined  from  displacement  measurements 
over  gage  lengths  that  included  a  number  of  mortar  joints.  Bedard  and 
Kotsovos  (1986)  presented  a  numerical  description  of  the  fracture  proc¬ 
esses  of  concrete  that  is  suitable  for  use  with  nonlinear  finite  element 
analysis  methods. 

Dhanasekar  and  Page  (1986)  studied  the  influence  of  brick  masonry  infill 
properties  on  the  behavior  of  infilled  frames,  using  a  finite  element  model 
to  simulate  the  behavior  of  infilled  frames  subjected  to  a  racking  load.  The 
finite  element  program  incorporated  a  material  model  for  the  infill  brick 
masonry  that  includes  appropriate  elastic  properties,  inelastic 
stress/strain  relations,  and  a  failure  surface.  Han  and  Chen  (1987)  sum¬ 
marized  the  efforts  in  developing  a  model  for  progressive  failure  analysis 
of  concrete  structures.  The  main  features  of  the  proposed  short-term, 
time-independent  constitutive  model  included  the  William  and  Warnke 
(1974)  five-parameter  or  Hsieh-Ting-Chen  (1988)  four-parameter  failure 
surface,  the  nonuniform  hardening  rule,  the  nonassociated  flow  rule  with 
changing  dilatancy  factor,  linear  softening  for  post-cracking  in  tension, 
and  multiaxial  softening  for  post-failure  in  compressions.  Rots  and  Borst 
(1987)  described  a  smeared  crack  model  that  covers  tensile  softening  in 
Mode  I  and  shear  softening  in  Mode  II  fracture.  In  addition,  the  model  ac¬ 
counted  for  unloading  and  reloading  and  for  multiple-crack  formation. 
Particular  forms  of  tension  and  shear  softening  functions  and  relations 
with  more  conventional  models  were  discussed.  Wiliam,  Pramono,  and 
Sture  (1987)  discussed  renewed  attempts  in  terms  of  fixed  and  rotating 


ERDC/CERL  TR-  08-19 


10 


crack  models  that  resort  to  fracture  mechanics  to  refine  the  traditional 
orthotropic  crack  formulation.  Following  a  similar  approach,  the  original 
concept  involving  Mode-I  type  cracking  is  being  broadened  to  include 
mixed-mode  fracture  interpretation  of  the  shear-retention  factor  if  the 
crack  memory  is  fully  retained.  Hsieh  et  al.  (1988)  discussed  an  elastic- 
plastic  fracture  constitutive  model  for  concrete  materials.  In  this  model,  a 
four-parameter  failure  criterion  was  used  to  define  the  ultimate  strength. 
The  associated  flow  rule  for  displaying  mixed  hardening  was  adopted  to 
represent  the  nonlinear  stress/strain  response.  A  crushing  coefficient  was 
introduced  to  identify  the  mode  of  failure,  and  an  anisotropic  elastic  rela¬ 
tionship  was  used  to  describe  the  post-failure  behavior.  Yazdani  and 
Schreyer  (1990)  developed  a  combined  plasticity  and  damage  mechanics 
model  for  concrete  within  the  general  framework  of  the  internal  variable 
theory  of  thermodynamics.  The  necessity  of  using  both  plasticity  and  dam¬ 
age  mechanics  was  discussed,  and  the  corresponding  surfaces  were  devel¬ 
oped  via  the  internal  dissipation  inequality. 

Lotfi  and  Shing  (1991)  assessed  the  capability  of  smeared  crack  models  in 
capturing  the  strength  and  various  failure  mechanisms  of  reinforced  ma¬ 
sonry  shear  walls.  A  smeared  finite  element  formulation  adopting  the  J2 
plasticity  model  for  uncracked  masonry  and  a  nonlinear  orthotropic  con¬ 
stitutive  model  for  crack  masonry  were  developed.  Stevens  and  Liu  (1992) 
combined  the  theories  of  continuum  damage  mechanics  and  plasticity  in  a 
strain-based  phenomenological  approach  to  yield  an  effective  constitutive 
model  for  plain  concrete.  Jamal  et  al.  (1992)  conducted  a  nonlinear  finite 
element  study  to  evaluate  the  in-plane  behavior  of  masonry  infilled  steel 
frames. 

Stankowski  et  al.  (1993)  presented  a  constitutive  theory  that  describes 
fracture  and  slip  of  an  interface  in  cementitious  composites  such  as  con¬ 
crete.  The  constitutive  model  was  developed  and  formulated  in  analogy  to 
the  theory  of  incremental  plasticity.  Adhesion,  debonding,  and  mobilized 
friction  mechanisms  were  considered  in  the  characterization  of  the  inter¬ 
face  behavior.  Abu-Lebdeh  and  Voyiadjis  (1993)  developed  a  model  that 
combined  plasticity  and  damage  mechanics  to  assess  both  multi-axial 
monotonic  and  cyclic  behavior  of  concrete.  The  model  adopted  a  bounding 
surface  concept  and  combined  plastic  deformation  with  the  deformation 
due  to  damage.  Lotfi  and  Shing  (1994)  simulated  the  failure  of  unrein¬ 
forced  masonry  panels  subjected  to  lateral  loads  by  means  of  a  finite  ele¬ 
ment  approach  in  which  the  mortar  joints  were  modeled  with  interface 
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elements  and  the  masonry  units  were  modeled  with  smeared  crack  ele¬ 
ments.  Lourenco  et  al.  (1995)  presented  two  models  for  the  micro  and 
macro  analysis  of  masonry  structures.  For  the  micro  modeling  of  masonry, 
an  interface  failure  criterion  that  included  a  straight  tension  cutoff,  the 
Coulomb  friction  law,  and  an  elliptical  cap  was  proposed.  Fuschi  et  al. 
(1995)  used  a  numerical  approach  for  the  structural  analysis  of  masonry 
walls  under  in-plane  stress  conditions.  Manzouri  et  al.  (1996)  presented 
two  elastic/viscoplastic  constitutive  models  for  the  modeling  of  masonry 
structures:  (1)  a  two-dimensional  continuum  model  based  on  the  Drucker- 
Prager  failure  criterion  for  modeling  masonry  units  and  (2)  an  interface 
model  for  modeling  mortar  joints.  Lee  et  al.  (1996)  introduced  a  homog¬ 
enization  technique  to  investigate  the  elastic-brittle  behavior  of  masonry 
panels  subject  to  incremental  lateral  loading.  Two  successive  steps  of  ho¬ 
mogenization  were  used  to  obtain  equivalent  elastic  properties  for  model¬ 
ing  the  elastic  behavior  of  masonry.  Barzegar  and  Maddipudi  (1997)  pre¬ 
sented  a  triaxial  constitutive  model  for  finite-element  analysis  of  pre-  and 
post-cracking  behavior  of  concrete.  A  five-parameter  ultimate  strength  en¬ 
velope  was  used  in  the  stress  space.  Lourenco  and  Rots  (1997)  evaluated 
the  performance  of  an  interface  elastoplastic  constitutive  model  for  the 
analysis  of  unreinforced  masonry  structures.  Both  masonry  components 
were  discretized,  aiming  for  a  rational  unit/joint  model  able  to  describe 
cracking,  slip,  and  crushing  of  the  material. 

Mehrabi  and  Shing  (1997)  carried  out  experimental  and  analytical  studies 
to  investigate  the  performance  of  masonry-infilled  R/C  frames  under  in¬ 
plane  lateral  loadings.  A  smeared-crack  finite  element  model  was  used  to 
model  the  behavior  of  concrete  in  the  R/C  frames  and  masonry  units.  It 
was  shown  that  finite  element  models  are  able  to  simulate  the  failure 
mechanisms  exhibited  by  infilled  frames,  including  the  crushing  and 
cracking  of  the  concrete  frames  and  masonry  panels,  and  the  sliding  and 
separation  of  the  mortar  joints.  Balan  et  al.  (1997)  presented  a  three- 
dimensional  concrete  material  model  for  nonlinear  finite-element  analysis 
of  concrete  solids  under  short-term  cyclic  loading.  The  proposed  model 
was  a  hypoplastic  orthotropic  model  based  on  a  stress-equivalent  uniaxial 
strain  relation  that  was  generalized  to  take  into  account  triaxial  stress 
conditions.  Carol  et  al.  (1997)  presented  a  simple  but  general  model  for 
normal  and  shear  cracking  in  quasi-brittle  materials.  The  model  was  de¬ 
fined  in  terms  of  the  normal  and  shear  stresses  on  the  average  plane  of  the 
crack,  and  the  corresponding  normal  and  shear  relative  displacements. 
Chiou  et  al.  (1999)  investigated  the  structural  behavior  of  a  framed  ma- 
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sonry  wall  subjected  to  in-plane  monotonic  loading  by  a  full-scale  test  and 
the  method  of  discontinuous  deformation  analysis.  The  concept  of  artifi¬ 
cial  joints  was  adopted  to  refine  discontinuous  deformation  analysis  so  it 
could  analyze  both  continuous  and  discontinuous  behavior  of  the  masonry 
structure.  Guinea  et  al.  (2000)  presented  a  micro-mechanical  model  for 
analyzing  Mode-I  fracture  of  brick  masonry.  The  analysis  was  based  on  a 
detailed  modeling  of  brick  and  mortar  fracture  by  means  of  the  fictitious 
(or  cohesive)  crack  model. 

Selection  of  constitutive  models 

The  behavior  of  reinforced  concrete  (R/C)  frames  infilled  with  unrein¬ 
forced  concrete  masonry  (UCM)  blocks  under  lateral  in-plane  loading  de¬ 
pends  strongly  on  the  development  of  diagonal  and  horizontal  cracks 
through  mortar  joints  and  the  evolution  of  these  cracks  under  shear  rever¬ 
sal,  and  repetitive  displacements.  Mehrabi  et  al.  (1994)  demonstrated 
through  experimental  and  analytical  studies  that  diagonal  and  horizontal 
cracking  within  the  infill  and  slip  at  cracked  joints  is  the  dominant  failure 
mechanism  for  R/C  frames  infilled  with  unreinforced  masonry.  Although 
cracking  and  crushing  and  accumulation  of  damage  also  occur  in  the  ma¬ 
sonry  units,  it  is  the  degradation  of  shear  resistance  of  cracked  masonry 
which  defines  degradation  of  lateral  stiffness  and  resistance  of  a  masonry 
infilled  frame.  Therefore,  if  a  constitutive  model  is  adopted  that  is  capable 
of  simulation  of  shear  degradation  along  masonry  joints  and  is  able  to 
model  quasi-brittle  material  such  as  concrete  and  masonry,  the  behavior 
of  an  infilled  frame  could  be  modeled  properly.  Hence,  the  model  selection 
was  based  on  capabilities  of  reviewed  models  in  simulation  of  masonry 
joint  behavior. 

Based  on  the  results  of  the  literature  review  reported  above,  it  is  concluded 
that  a  combination  of  continuum  and  interface  constitutive  models  are  the 
most  appropriate  for  R/C  frames  infilled  with  UCM.  The  continuum  model 
is  used  for  R/C  in  frame  and  masonry  units  in  infill,  and  the  interface 
model  is  used  for  mortar  joints  between  masonry  units  as  well  as  joints  be¬ 
tween  the  infill  and  frame. 

Characteristics  and  application  of  selected  constitutive  models 

To  model  the  concrete  in  the  R/C  frames  and  the  masonry  units  in  the  in¬ 
fill  panels,  the  smeared-crack  finite  element  formulation  developed  by 
Lotfi  and  Shing  (1991)  was  selected.  In  this  formulation,  a  J2-plasticity 
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model  with  an  isotropic  strain-hardening/softening  law  was  utilized  to 
model  an  uncracked  material.  The  plasticity  model  is  combined  with  the 
Rankine  tension-cutoff  criterion  to  signal  the  onset  of  cracking. 

To  model  the  behavior  of  the  mortar  joints  between  masonry  units  as  well 
as  the  joints  at  the  infill-to-frame  interfaces,  the  interface  model  developed 
by  Mehrabi  and  Shing  (1995)  was  selected. 

The  longitudinal  reinforcing  bars  in  the  R/C  frame  members  can  be  mod¬ 
eled  as  one-dimensional  (lD)  truss  elements  using  an  elastic-hardening- 
plastic  material  model  available  in  all  commercial  FE  programs.  The  shear 
reinforcement  in  R/C  frame  members  can  be  modeled  as  two-dimensional 
(2D)  shell  elements  using  elastic-hardening-plastic  material  model.  The 
reinforcement  can  be  connected  to  concrete  elements  assuming  a  perfect 
bond,  or  through  bond-slip  FE  elements  such  as  those  developed  by  Me¬ 
hrabi  and  Shing  (1995).  However,  Mehrabi  and  Shing  demonstrated  that 
the  use  of  bond-slip  elements  between  concrete  and  reinforcement  at 
nodal  points  has  a  minimal  effect  on  modeling  the  behavior  of  infilled 
frames.  Therefore,  the  use  of  bond-slip  elements  for  reinforcement  does 
not  seem  necessary. 

The  smeared-crack  model 

The  smeared-crack  model  developed  by  Lotfi  and  Shing  (1991)  uses  a  J2  - 
plasticity  model  with  an  isotropic  strain-hardening/softening  law  for  mod¬ 
eling  an  uncracked  material.  The  plasticity  model  is  combined  with  the 
Rankine  tension-cutoff  criterion  to  signal  the  onset  of  cracking.  The  failure 
criteria  are  shown  in  terms  of  the  principal  stresses,  ai  and  02,  in  Figure  2. 
The  isotropic  strain-hardening/softening  is  governed  by  ae  =  <re  (s  ) , 

which  consists  of  a  parabolic  curve  and  an  exponential  tail  as  shown  in 
Figure  2d.  A  crack  initiates  when  the  maximum  principal  stress  reaches 
the  tensile  strength,  f  't ,  in  a  direction  normal  to  the  maximum  principle 

stress.  The  cracked  material  is  assumed  to  be  nonlinear  orthotropic,  with 
the  axes  of  orthotropy,  n-t,  normal  and  tangential  to  the  direction  of  the 
crack  (Figure  2b).  For  a  coaxial  rotating  crack  model,  which  is  considered 
for  this  model,  the  crack  direction  and  the  axes  of  orthotropy  rotate  with 
the  principal  axes  of  strain  in  such  a  way  that  the  crack  remains  normal  to 
the  direction  of  the  maximum  principal  strain.  Variations  of  tension  stress 
and  tension  stiffness  before  and  after  cracking  are  shown  in  Figure  2c. 
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(a)  Failure  Criteria 


(b)  Crack  Orientation 


e!P  e2P 


(d) 


Effective  Stress-Effective 
Plastic  Strain  Curve 


(e) 


E,  E2  El 

Uniaxial  Behavior  Parallel 
to  the  Crack 


Figure  2.  Smeared  crack  model  failure  criteria. 


For  the  case  of  a  coaxial  rotating  crack,  the  shear  modulus  is  given  by 

Gnt  =  >  where  cti  and  a2  are  the  principal  stresses,  and  si  and  so 

2(£l-£2) 

are  the  principal  strains  (Bazant  1983).  The  tangent  modulus  in  the  direc¬ 
tion  parallel  to  the  crack  is  represented  by  Ett.  The  compression  behavior 
in  a  direction  parallel  to  the  crack  is  modeled  by  a  nonlinear  harden¬ 
ing/softening  law  similar  to  the  uniaxial  compressive  behavior  given  by 
the  plasticity  model  for  an  uncracked  material.  The  uniaxial  effective  stress 
versus  effective  plastic  strain  and  uniaxial  behavior  parallel  to  the  crack 
are  shown  in  Figure  2d  and  9e,  respectively.  Additional  details  of  the  crack 
model  can  be  found  in  Lotfi  and  Shing  (1991). 

The  dilatant  interface  model 

To  model  the  behavior  of  concrete,  rock,  and  mortar  joints,  various  plastic¬ 
ity-based  interface  constitutive  models  have  been  proposed  by  Plesha  et  al. 
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(1989),  Stankowski  (1990),  Stankowski  et  al.  (1993),  Carol  and  Prat  (1991), 
and  Lotfi  and  Shing  (1994).  However,  there  are  still  some  aspects  of  the 
physical  behavior  of  interfaces  that  have  not  been  considered  properly  in 
most  existing  models.  These  include  the  compressive  hardening  behavior 
of  interfaces,  the  reversal  of  shear  dilatancy  in  the  case  of  cyclic  loading, 
and  the  normal  contraction  of  an  interface  under  shear  sliding.  These 
physical  properties  of  interfaces  are  incorporated  in  the  constitutive  model 
proposed  by  Mehrabi  and  Shing  (1997). 

In  this  model,  the  nonlinear  hardening  behavior  of  an  interface  under  a 
compressive  stress  is  assumed  to  be  an  elastic  phenomenon.  No  dilatation 
is  assumed  in  the  elastic  regime.  In  the  plastic  regime,  the  total  shear  dila¬ 
tation  exhibited  by  an  interface  is  considered  to  be  the  sum  of  a  normal 
compaction  and  geometric  dilatation.  The  latter  is  due  to  the  wedging  ac¬ 
tion  of  the  asperities,  as  shown  in  Figure  3.  Owing  to  the  reversible  nature 
of  this  dilatation,  the  dilatation  is  conceived  as  a  geometric  phenomenon. 
The  normal  contraction  of  an  interface,  due  to  the  loss  of  loose  particles,  is 
an  irreversible  phenomenon  that  is  modeled  using  plasticity  theory. 


A  three-parameter  hyperbolic  yield  criterion  proposed  by  Lotfi  and  Shing 
(1994)  for  cohesive  interfaces  is  adopted  here,  as  shown  in  Figure  4.  As  il¬ 
lustrated,  evolution  and  shrinkage  of  the  yield  surface  occurs  with  the  oc¬ 
currence  of  plastic  work. 
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Figure  4.  Hyperbolic  yield  criterion  for  cohesive  interfaces. 


To  model  the  normal  compaction  of  an  interface  due  to  the  loss  of  dam¬ 
aged  materials  as  well  as  the  opening  of  a  fractured  interface,  a  non- 
associated  flow  rule  is  proposed,  as  shown  in  Figure  5.  The  direction  of  the 
plastic  relative  displacement  is  governed  by  the  flow  rule. 
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Figure  5.  Plastic  potential  (non  associated  flow  rule). 


The  shear  dilatancy  is  considered  as  the  sum  of  geometric  dilatation  and 
normal  plastic  compaction.  Geometric  dilatation  occurs  only  after  fracture, 
and  it  is  expressed  as  a  function  of  the  plastic  shear.  The  interface  is  as¬ 
sumed  to  function  as  an  asperity,  the  slope  of  which  decreases  with  respect 
to  the  cumulative  plastic  work  performed  by  shear.  The  loading  and 
unloading  satisfy  the  Kuhn-Tucker  conditions.  For  the  integration  of  the 
above  constitutive  relations,  the  return  mapping  algorithm  based  on  the 
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operator-splitting  method  proposed  by  Ortiz  (1981)  is  adopted.  The  solu¬ 
tion  of  the  rate  equations  is  decomposed  into  an  elastic  prediction  and  a 
plastic  correction.  In  the  elastic  prediction,  the  plastic  strains  and  internal 
variables  are  assumed  to  remain  constant,  and  the  applied  displacement 
increment  is  assumed  to  be  totally  elastic.  For  the  plastic  correction,  there 
is  no  additional  displacement  increment. 


ERDC/CERL  TR-  08-19 


18 


3  Software  Implementation  of  Constitutive 
Models 

Constitutive  model  performance  in  previous  FE  analysis 

Mehrabi  et  al.  (1994)  and  Mehrabi  and  Shing  (1997)  examined  the  validity 
of  the  constitutive  models  by  using  them  in  a  finite  element  analysis  based 
on  previous  experimental  testing  of  physical  masonry-infilled  R/C  frame 
models.  Specimens  with  weak  and  strong  frame  designs  were  analyzed  for 
performance  when  subjected  to  a  constant  vertical  compressive  load.  The 
concrete  frames  were  modeled  with  nine-node  quadrilateral  smeared  crack 
elements,  and  the  shear  reinforcement  was  smeared  out  over  concrete 
elements.  The  longitudinal  bars  in  the  frame  were  modeled  with  two-node 
bar  elements.  They  were  connected  to  the  nine-node  concrete  elements  at 
the  two  external  nodes.  For  masonry  units,  four-node  smeared  crack  ele¬ 
ments  were  used.  The  mortar  joints  in  the  masonry  panels  and  along  the 
interfaces  between  an  infill  and  the  frame  were  modeled  by  two  double¬ 
node  interface  elements.  To  model  the  possible  shear  failure  of  the  col¬ 
umns,  three  double-node  interface  elements  were  used  at  critical  locations 
near  the  top  and  bottom  sections  of  the  columns.  Extensive  materials  test¬ 
ing  was  also  carried  out.  The  constitutive  models  were  calibrated  with  the 
results  of  the  physical  modeling  reported  in  Mehrabi  et  al.  (1994). 

FE  software  selection 

It  was  decided  to  use  a  commercial  FE  software  program  in  this  study  for 
purposes  of  economy,  ease  of  model  construction,  inclusion  of  user- 
friendly  input  and  output  formats,  and  integrated  graphics  capabilities. 
Four  finite  element  packages,  each  with  capabilities  for  modeling  concrete 
elements  and  interfaces,  were  reviewed.  They  were  ANSYS,  ADINA, 
ABAQUS,  and  DIANA.  The  programs  were  reviewed  to  understand  their 
capabilities  in  modeling  structural  discontinuities  such  as  mortar  joints 
and  modeling  a  concrete  and  masonry  continuum.  For  the  type  of  infilled 
frames  considered  in  this  study,  the  defining  parameters  are  cracking  and 
separation  of  mortar  joints  and  infill -to-frame  interfaces,  and  cracking  and 
crushing  of  the  infill  material.  Degradation  of  shear  properties,  especially 
at  the  interfaces  and  joints,  are  of  specific  interest.  Based  on  those  criteria, 
DIANA  identified  as  the  best  fit  for  the  purposes  of  this  study. 
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Capabilities  of  the  DIANA  program 

This  section  examines  the  capabilities  and  limitations  of  DIANA  in  terms 
of  analyzing  masonry-infilled  R/C  frames  subjected  to  in-plane  lateral 
loading. 

General  modeling 

In  DIANA,  the  analytical  model  of  a  structure  can  be  divided  into  three 
major  parts:  (l)  the  discretization  of  the  geometry  through  a  finite  element 
mesh,  (2)  the  modeling  of  the  physical  behavior  of  the  materials  used  in 
the  structure,  such  as  concrete  and  reinforcing  steel,  and  (3)  the  modeling 
of  the  structural  effects  that  influence  the  behavior  of  the  structure,  such 
as  large  displacements  or  time-dependent  behavior. 

Structural  behavior 

To  make  the  discretization  in  DIANA,  the  user  needs  knowledge  of  applied 
mechanics  and  finite  elements  in  particular,  but  it  is  assumed  for  this  part 
that  the  geometry  of  the  structure  has  been  discretized  in  a  proper  fashion. 
For  the  modeling  of  the  physical  behavior  of  the  material,  the  user  needs 
knowledge  about  failure  mechanisms  that  can  occur  in  the  structure.  In 
concrete  structures,  for  example,  the  behavior  is  mainly  influenced  by 
cracking  and  crushing  of  the  concrete  and  yielding  of  the  reinforcement. 
Structural  effects  are  those  phenomena  that  can  influence  the  behavior  of 
the  structure,  and  they  depend  largely  on  the  history  of  the  structure. 

Modeling  brittle  materials 

To  model  concrete  structures,  or  structures  made  of  brittle  and  quasi- 
brittle  materials  in  general,  DIANA  offers  a  broad  range  of  element  types. 
The  constitutive  behavior  of  quasi-brittle  material  is  characterized  by  ten¬ 
sile  cracking  and  compressive  crushing,  and  by  long-term  effects  like 
shrinkage  and  creep. 

Cracking 

Cracking  can  be  modeled  in  DIANA  with  a  multidirectional  fixed-crack 
model  with  tension  softening  and  shear  retention.  Brittle  cracking,  linear 
tension  softening,  multilinear  softening,  and  nonlinear  softening  accord¬ 
ing  to  Hordijk  (1991)  is  available.  Also  a  plasticity-based  formulation  for 
cracking  is  available:  the  principal  Rankine  stress  criterion  resembles  the 
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rotating  crack  model.  However,  this  model  is  applicable  only  to  plane 
stress,  plane  strain,  and  axisymmetric  elements. 

In  multi-axial  stress  states,  the  compressive  stress  accounted  for  in  the 
model  can  exceed  the  compressive  strength  of  the  material.  In  such  a  case, 
the  crack  model  can  be  combined  with  a  plasticity  model  that  describes  the 
crushing  of  the  material.  The  Mohr-Coulomb  and  Drucker-Prager  models 
are  especially  applicable  to  quasi-brittle  structures. 

The  combination  of  tensile  and  compressive  stresses  can  also  be  modeled 
in  DIANA  with  a  multi-surface  plasticity  model  that  is  available  for  biaxial 
stress  states.  However,  this  model  too  is  applicable  only  to  plane  stress, 
plane  strain,  and  axisymmetric  elements. 

Reinforcement 

The  reinforcement  in  a  concrete  structure  can  be  modeled  with  the  em¬ 
bedded  reinforcement  types  available  in  DIANA.  The  constitutive  behavior 
of  the  reinforcement  can  be  modeled  by  an  elastoplastic  material  model 
with  hardening.  Furthermore,  temperature  influence  on  the  material’s 
Young's  modulus  and  the  thermal  expansion  coefficient  can  be  taken  into 
account. 

Specific  support  for  modeling  masonry 

Masonry  structures  are  analyzed  on  two  different  levels:  the  macro  level, 
where  global  behavior  is  simulated;  and  the  meso  level,  where  behavior  is 
analyzed  in  more  detail.  For  macro-level  analysis,  DIANA  offers  multi¬ 
directional  fixed-crack  models  and  plasticity  models  to  simulate  cracking 
and  crushing,  respectively.  However,  the  orthotropic  nature  of  masonry 
cannot  be  modeled  with  the  fixed  crack  and  standard  plasticity  models  be¬ 
cause  these  models  involve  isotropic  elasticity,  and  they  do  not  allow  com¬ 
bination  with  orthotropic  elasticity.  The  anisotropic  Rankine-Hill  plastic¬ 
ity  model  is  appropriate  for  modeling  masonry  because  it  allows 
orthotropic  elasticity  to  be  employed  and  incorporates  different  strength 
and  degradation  parameters  to  simulate  the  different  behaviors  parallel 
and  perpendicular  to  the  masonry  bed. 

For  meso-level  analysis,  the  orthotropy  of  masonry  blocks  can  be  modeled 
by  continuum  elements,  and  the  joints  can  be  modeled  by  interface  ele¬ 
ments.  For  this  type  of  task  there  are  various  models  to  describe  the  inter- 
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face  behavior,  including  a  discrete  crack  model,  a  Coulomb  friction  model, 
and  a  combined  Coulomb  friction/tension  cutoff/compression  cap  model. 

The  models  described  above  are  not  the  only  models  that  can  be  applied  to 
masonry.  In  some  cases,  inclusion  of  the  elastic  orthotropy  of  the  masonry 
may  not  be  essential,  and  a  standard  smeared  or  total  strain  crack  model 
with  isotropic  elasticity  may  be  applied  as  well.  That,  for  instance,  is  the 
case  when  the  orthotropy  factor  is  low  or  when  cracks  are  expected  to  oc¬ 
cur  predominantly  in  only  one  direction.  Also,  the  combined  fric¬ 
tion/tension/compression  interface  model  is  not  always  required,  and  one 
may  choose  to  use  a  standard  discrete  crack  or  Coulomb  friction  model.  It 
should  be  realized  that  the  Rankine-Hill  continuum  model  and  the  com¬ 
posite  interface  model  are  plasticity  based,  which  means  that  unloading 
takes  place  elastically.  To  model  secant  unloading  and  crack  clo¬ 
sure/reopening  specifically,  other  constitutive  models  may  be  chosen. 

Verification  of  constitutive  models  in  DIANA 

Procedure 

For  the  purpose  of  this  study,  it  was  decided  to  perform  modeling  and 
analysis  at  the  meso  level  as  follows:  (l)  masonry  units  and  R/C  frames 
were  modeled  using  continuum  elements  with  smeared-crack-type  consti¬ 
tutive  models  capable  of  cracking  and  crushing,  and  (2)  mortar  joints  were 
modeled  using  structural  interfaces  capable  of  modeling  cohesion,  separa¬ 
tion,  shear  degradation,  cyclic  behavior,  and  closing. 

In  order  to  evaluate  how  the  finite  element  models  of  interest  in  this  inves¬ 
tigation  perform  within  DIANA,  it  was  necessary  to  observe  the  behavior 
of  those  models  under  different  modes  of  deformation.  The  selected  inter¬ 
face  analytical  model  for  the  mortar  joints  was  calibrated  by  evaluating  it 
against  results  from  an  established  discrete  crack  interface  model  (here 
called  the  comparative  model).  Displacement  comparisons  between  the 
two  models  were  made  normal  to  the  interface,  parallel  to  the  interface, 
and  in  combinations  of  the  two.  The  solid  masonry  material  was  modeled 
using  a  smeared  crack  model  in  order  to  demonstrate  different  modes  of 
progressive  failure.  The  basic  interface  and  solid  smeared  crack  definitions 
and  implementations  were  shown  previously,  in  Figure  2. 

Prior  to  performing  evaluations  and  comparisons,  the  intended  capabili¬ 
ties  of  each  analytical  model  were  listed  in  a  table  and  compared.  Table  1 
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shows  those  basic  comparisons  for  the  interface  model.  The  interface 
model  built  into  DIANA  includes  a  compression  cap  that  can  represent 
compressive  damage  associated  with  the  failure  in  joints,  but  the  compara¬ 
tive  model  does  not  have  this  feature.  That  difference  may  show  up  in  a 
three-dimensional  analysis.  Another  difference  is  that  the  DIANA  interface 
model  does  not  represent  cumulative  damage  to  the  joints  in  the  form  of 
loss  of  material  while  the  comparative  model  does  have  this  feature.  This 
cumulative  damage  feature  affects  the  active  dilatancy,  and  it  can  be  im¬ 
portant  once  the  joints  are  placed  under  constraints  such  as  an  infill 
frame.  (A  future  study  will  address  that  modeling  capability  in  more  de¬ 
tail.) 

Despite  the  two  significant  differences,  the  two  interface  models  share  a 
wide  basis  for  representing  frictional  behavior,  frictional  and  tensile  deg¬ 
radation,  and  progress  of  damage  to  the  properties,  and  they  may  be  ex¬ 
pected  to  behave  similarly.  Features  and  parameters  that  both  models 
share  are  shown  in  Table  l. 


Table  1.  Comparison  of  DIANA  model  and  comparative  model. 


Model  Setup 
Description 

Shape 

DIANA  model 

Comparative  model 

Comment 

Appendix 

Equation 

Parameters 

Appendix 

Equation 

Parameters 

Shear  Slipping 

Coulomb 

A1 

4 

©,  c 

A2 

5 

|d,s,r 

similar 

4a 

c 

13 

p 

p=<5 

5 

© 

13 

s 

13 

r 

r=(c2-|u.2s2)/2s 

Flow  Rule 

Plastic  Potential 

A1 

13,14 

Y,  5,au 

A2 

6 

a,n 

Tension 

softening 

Exponential 

A1 

16 

Y,5 

A2 

13 

S 

8,9 

13 

r 

Tension  Cut-Off 

Cut-Off  Line(DIANA) 

A1 

17,18 

ft 

- 

- 

- 

Tip  of  Hyperbola 

A2 

5 

s,r 

Compression 

Cap 

Elliptic 

A1 

22 

Co 

NA 

NA 

NA 

au 

Dilatancy 

- 

A1 

12 

Y,  5,au 

A2 

8,12 

S.Y 

Loss  of  Material 

- 

NA 

NA 

NA 

A2 

9 

5.Y 

Verification  and  calibration  of  interface  model 

For  verification,  the  interface  model  was  tested  for  its  overall  behavior.  A 
single  joint  placed  in  between  two  elastic  masonry  blocks  was  tested  under 
compression,  tension,  and  direct  shear  loads.  The  methods  of  analysis, 
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namely  Modified  Newton-Raphson  and  the  secant  method  with  BFGS 
(Broyden,  Fletcher,  Goldfarb,  and  Shannoare),  are  used  in  different  analy¬ 
ses  and  checked  for  efficiency.  The  direct  solver  is  employed,  and  the  con¬ 
vergence  criterion  is  set  for  a  small  value  of  the  norm  of  the  residual  en¬ 
ergy  at  each  iteration.  The  step  size  is  varied  for  each  run  in  order  to  find 
and  fit  the  parameters  to  their  best  performance. 

Verification  and  calibration  analyses 

For  verification  purposes,  three  different  analyses  of  an  interface  element 
were  performed  using  DIANA.  The  first  analysis  was  executed  to  check  the 
tensile  capacity  of  the  interface  element,  as  shown  in  Figure  6.  The  inter¬ 
face  element  reached  a  tensile  capacity  of  38.2  psi  before  significant  joint 
opening  and  tension  softening  occurred.  The  second  analysis  consisted  of 
a  shear  softening  test  with  no  applied  normal  pressure,  as  represented  by 
Figure  7.  The  third  analysis  was  a  shear  softening  test  with  an  applied 
normal  pressure  of  100  psi,  as  shown  in  Figure  8.  As  expected,  the  normal 
compressive  force  increased  the  shear  capacity,  yielding  a  maximum  shear 
stress  of  123.8  psi.  Figure  7  and  Figure  8  demonstrate  the  DIANA  imple¬ 
mented  constitutive  model.  The  input  values  for  the  interface  element  ma¬ 
terial  parameters  are  given  in  Table  2,  which  is  followed  by  a  definition  of 
the  parameter  symbols.  The  behavior  of  the  joint  under  a  constant  normal 
pressure  of  100  psi  followed  by  a  cyclic  shear  load  was  tested  for  one  full 
cycle,  as  shown  in  Figure  9.  For  demonstration  purposes,  Figure  10  is 
shown  to  represent  the  same  plot  as  Figure  9  using  the  DIANA  post¬ 
processor. 


Table  2.  Material  parameters  for  DIANA  interface  element  verification  analysis. 


Knn 

Kss 

ft 

G/ 

c 

<|>i 

Yi 

<K 

CTu 

5 

Gt" 

f'c 

Cs 

Gfc 

Kp 

2.8E+5 

3.5E+5 

40 

1.69 

40 

0.79 

0.005 

0.65 

150 

2.3 

16.1 

1,500 

1.0 

2 

0.006 

Definition  of  symbols: 

Knn  and  Kss  are  the  normal  and  tangential  numerical  stiffness  parameters  of  the  joint 

ft  =  joint  tensile  strength 

Gt1  =  first  mode  fracture  energy 

c  =  cohesion 

(j)i  =  internal  friction  angle 

\|/  =  dilatancy  angle 

tp0  =  initial  internal  friction  angle 

4>r  =  residual  internal  friction  angle 

ctu  =  confining  stress  over  which  the  dilatancy  will  be  zero 

5  =  dilatancy  shear  slip  degradation  coefficient 
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Gf11  =  shear  mode  fracture  energy 
f'c  =  compressive  strength 

Cs  =  parameter  controlling  the  shear  stress  contribution  to  failure 
Gfc=  compressive  failure  fracture  energy 
kp  =  norm  of  plastic  strain 


Joint  opening  (in.) 


Figure  6.  Interface  element  tension  softening  curve  from  DIANA  analysis. 


Shear  Displacement  (in.) 


Figure  7.  Interface  shear  softening  under  zero  normal  stress  from  DIANA  analysis. 
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Shear  Displacement  (in.) 

8.  Interface  shear  softening  under  100  psi  normal  stress  from  DIANA  analysis. 
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Figure  9.  Interface  element  subjected  to  cyclic  shear  loading  under  100  psi  normal  stress. 
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i DIANA  8.1-460  :  U.S.Army  Construction  Eng. Res.  22-MAY-2004  23:03  trial  1 


MAXXMIN  ON  WHOLE  GRAPH: 
YMAX  =  297  YMIN  =  -445 
XMAX  =  1.2  XMIN  =  -1.2 
VARIATION  OVER  LOADCASES 
NODE  6 


Figure  10.  Interface  element  subjected  to  cyclic  shear  loading  under  100  psi  normal  stress, 

as  output  by  the  DIANA  post-processor. 

Numerical  verification  of  the  interface  model 

The  capability  of  the  interface  constitutive  model  was  validated  using  data 
produced  in  the  direct  shear  tests  conducted  on  mortar  joints  performed 
by  Mehrabi  et  al.  (1994).  The  constitutive  model  was  implemented  in  a 
one-dimensional  isoparametric  interface  element.  For  the  calibration  of 
the  model,  data  from  shear  tests  conducted  under  three  different  levels  of 
normal  stresses  (100  and  150  psi)  and  subjected  to  cyclic  loading  are  de¬ 
sirable. 

The  model  was  calibrated  based  on  the  parameters  of  the  comparative  in¬ 
terface  model.  These  characteristics  are  very  similar  to  parameters  defined 
previously  in  Table  2.  The  dilatancy  in  both  models  does  well  in  the  mono¬ 
tonic  loading  regime.  As  the  loads  are  changed  to  cyclic,  the  cumulative 
damage  represented  by  the  loss  of  joint  material  is  not  modeled  appropri¬ 
ately  by  the  standard  interface  model  built  into  DIANA.  The  loss  is  a  physi¬ 
cal  fact  observed  in  direct  shear  tests  on  masonry  bed  joints  and  rock 
joints. 


To  calibrate  and  verily  the  interface  analytical  model,  two  trials  were  run 
in  DIANA  with  a  slight  variation  in  the  material  parameters  for  the  inter- 
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face  element,  namely  5.  Table  3  details  the  material  parameters  of  the 
comparative  model  and  Table  4  details  the  parameters  used  for  the  two 
different  DIANA  trials. 

Next,  the  DIANA  models  were  run,  and  the  results  were  compared  with 
the  numerical  and  experimental  results  of  the  comparative  model  using 
the  two  different  trials  previously  noted.  Figure  11  represents  a  sample 
shear  deformation  of  the  interface  element  as  shown  in  the  DIANA  post¬ 
processing  window.  The  shear  stress  versus  shear  displacement  results  of 
the  calibrated  DIANA  model,  subjected  to  a  too  psi  normal  stress,  are 
shown  in  Figure  12.  Those  results  are  verified  using  the  results  of  the  com¬ 
parative  model  under  the  same  loading  conditions,  as  shown  in  Figure  13. 
The  results  demonstrate  the  model’s  capability  to  simulate  the  cyclic  be¬ 
havior  of  rock  and  mortar  joints  under  a  constant  normal  pressure. 


Table  3.  Material  parameters  of  the  comparative  model. 


Knn  (ksi) 

Dtt  (ksi) 

So  (psi) 

G/  (lb-in) 

c  (psi) 

<l>i 

<1)1- 

a 

5 

G^'  (lb-in) 

P 

ro  (psi) 

r<  (psi) 

So 

Y 

h 

28 

35 

40 

1.61 

26 

0.9 

0.8 

4 

0 

16.1 

2 

10 

5 

0.45 

3E-04 

3 

55 

Table  4.lnterface  element  material  parameters  for  trials  1  and  2. 


Trial 

Knn 

Kss 

ft 

Gf' 

c 

<t>i 

Ti 

<t>r 

CJu 

5 

Gt" 

f'c 

Cs 

Gtc 

Kp 

1 

2.8E+4 

3.5E+4 

40 

1.69 

26 

0.9 

0.005 

0.8 

150 

2.3 

16.1 

1,500 

1.0 

2 

0.006 

2 

2.8E+4 

3.5E+4 

40 

1.69 

26 

0.9 

0.005 

0.8 

150 

0.5 

16.1 

1,500 

1.0 

2 

0.006 
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iDIANA  8.1-460  :  U.S.Arny  Construction  Eng. Res. 


23-MAY-2Q04  20:52  directShear 


model:  diana_w 

LC2:  LOAD  CASE  2 
STEP:  25  LOAD:  .75E-1 
NODAL  TOTX...G  TDTX 
MAX  =  .75E-1  MIN  =  0 
FACTOR  =  10.6 


Figure  11.  Shear  deformation  of  the  interface  element  shown  as  DIANA  output. 


Figure  12.  Calibrated  DIANA  interface  element  (100  psi  normal  stress, au=150  psi). 
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Shear  displacement  (cm) 

Figure  13.  Comparative  model  results  for  shear  under  a  100  psi  normal  stress. 

The  results  from  DIANA  match  the  results  from  the  comparative  model 
very  closely.  In  Figure  14,  the  dilatancy  of  the  interface  element  under  a 
100  psi  normal  stress  is  compared  with  the  results  of  the  comparative 
model.  Model  results  for  dilatancy  under  100  psi  normal  stress  are  shown 
in  Figure  15. 
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Figure  14.  DIANA  window  showing  dilatancy  under  100  psi  normal  stress  (au=150  psi). 
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Figure  15.  Comparative  model  results  for  dilatancy  under  100  psi  normal  stress. 


The  calibrated  DIANA  interface  element  under  150  psi  normal  stress  is 
shown  in  Figure  16,  and  the  corresponding  dilatancy  under  150  psi  normal 
stress  is  shown  in  Figure  17.  Comparative  model  results  for  dilatancy  un¬ 
der  150  psi  normal  stress  are  shown  in  Figure  18. 


Shear  Displacement  (in) 

Figure  16.  Calibrated  DIANA  interface  element  under  150  psi  normal  stress. 


Vertical  Displacement  (in) 
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Verification  and  calibration  of  smeared  crack  model 

In  this  section,  the  primary  focus  is  on  the  verification  of  the  DIANA  total 
strain  smeared  crack  model  for  concrete  or  clay  masonry  blocks.  Without 
loss  of  generality,  a  plain  stress  condition  is  considered  in  this  evaluation 
to  match  the  two-dimensional  interface  model  described  previously  be¬ 
cause  the  two  models  will  be  used  together  to  model  masonry  assem¬ 
blages. 

To  verily  the  overall  behavior  of  the  smeared  crack  model,  a  single  block 
was  tested  under  compression,  tension,  and  direct  shear.  A  rotating 
smeared  crack  was  used,  for  which  the  retention  factor  G=i.o.  The  mate¬ 
rial  parameters  for  the  CMU  element  tested  under  these  conditions  are 
shown  in  Table  5,  followed  by  definitions  of  the  symbols.  The  methods  of 
analysis  and  the  convergence  rate  are  varied  for  each  run  in  order  to  find 
and  fit  the  parameters  to  their  best  performance. 


Table  5.  Material  parameters  for  a  single-unit  CMU  element  test. 


E 

V 

ft 

GH 

Tensile 

Curve 

Shear* 

Curve 

P 

f'c 

Comp. 

Curve 

Gc 

Gr" 

fc 

Cs 

GfC 

Kp 

3.5E+6 

0.16 

200 

0.03 

Exponential 

BEDIAG 

0.02 

3,000 

Parabolic 

22 

16.1 

1500 

1.0 

2 

0.006 

*  Shear  is  not  defined  for  a  rotating  crack  model. 


Definition  of  symbols: 

E  =  modulus  of  elasticity,  psi 

v  =  Poisson’s  ratio 

ft  =  tensile  strength,  psi 

Gt1  =  first  mode  fracture  energy 

P  =  the  shear  retention  coefficient 

f'c  =  the  compressive  strength,  psi 

Gc  =  the  fracture  energy  in  compression 

REDCRV  VC1993  =  Effect  of  lateral  cracking  on  compressive  strength 
CNFCRV  VECCHI  =  Effect  of  lateral  confinement  on  compressive  strength 

Compressive  behavior 

A  set  of  parameter  study  cases  was  run  on  a  four-noded  element  repre¬ 
senting  a  concrete  block.  In  this  study,  we  were  focused  on  the  compres¬ 
sive  fracture  energy  and  its  effect  on  the  parabolic  curve  of  compressive 
failure.  The  concrete  masonry  unit  (CMU)  model  was  calibrated  to  a  value 
of  Gf1  =  55,  f  m  =  3,000  psi,  and  a  modulus  of  elasticity,  E,  of  2,000,000  psi. 
The  nonlinear  behavior  of  the  CMU  element  subjected  to  axial  compres¬ 
sion,  as  shown  in  Figure  19,  exhibits  a  parabolic  stress/strain  curve. 
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Figure  19.  Compressive  elastic/plastic  hardening/softening  for  a  CMU. 


Tensile  behavior 

The  model  also  was  studied  for  an  uniaxial  tension  case.  No  reference  ex¬ 
periment  was  used  for  comparison  in  this  case,  so  the  model  was  just  cali¬ 
brated  for  reasonable  values  of  tensile  strain  associated  with  the  peak 
stress  (200/2E6=0.000i).  The  model  behaved  consistently,  and  the  algo¬ 
rithms  appeared  to  be  very  robust.  Figure  20  shows  pre-  and  post-peak 
results  from  a  monotonic  tensile  loading  of  a  CMU  block,  which  indicate 
an  exponential  curve  for  the  tensile  behavior.  The  model  was  calibrated  in 
tension  to  a  tensile  stress  of  200  psi  and  Gf1  =  18. 


Figure  20.Tensile  elastic/plastic  softening  of  a  CMU  element. 
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4  Modeling  and  Analysis  of  Masonry  Prisms 

Physical  modeling  of  Mehrabi  masonry  prisms 

The  three-unit,  one-half-scale  masonry  prisms  tested  by  Mehrabi  et  al. 
(1994)  consisted  of  concrete  hollow  or  solid  blocks  with  nominal  dimen¬ 
sions  of  4  in.  by  4  in.  by  8  in.  (3.625  in.  by  3.625  in.  by  7.625  in.  actual) 
with  3/8  in.,  Type  S  masonry  mortar  joints.  The  mortar  was  applied  onto 
face  shells  only  for  the  case  of  prisms  built  with  hollow  blocks,  and  onto 
the  entire  bed  joints  for  the  case  of  solid  block  prisms.  Hollow  blocks  had  a 
face  shell  thickness  of  0.625  in.  Therefore,  the  equivalent  width  of  mortar 
joint  for  hollow  blocks  was  assumed  to  be  twice  the  shell  thickness,  i.e., 
1.25  in.  An  equivalent  thickness  of  1.8  in.  for  hollow  blocks,  including  face 
shells  and  webs,  was  calculated  based  on  proportion  of  net  concrete  cross- 
section  with  respect  to  the  gross  cross-sections  of  the  block.  Table  6  sum¬ 
marizes  the  average  material  properties  and  test  results  for  hollow  and 
solid  masonry  prisms.  Calculations  are  based  on  the  results  reported  by 
Mehrabi  et  al.  (1994). 


Table  6.  Average  material  properties  and  axial  test  results 
for  masonry  prisms  (Mehrabi  et  al.  1994). 


Type  of  Blocks 

feu  (psi) 

f  cm  (psi) 

f  m  (psi) 

su  (in./in.) 

Hollow 

2,400 

2,100 

1,550 

0.0031 

Solid 

2,300 

2,100 

1,930 

0.0027 

This  table  includes  the  average  compressive  strength  of  masonry  units 
with  respect  to  the  net  cross-sectional  area  (feu),  average  compressive 
strength  of  masonry  mortar  tested  on  cubes  (f  cm),  compressive  strength  of 
masonry  prisms  with  respect  to  the  net  cross-sectional  area  (f  m),  and 
strain  at  maximum  strength  for  masonry  prisms  (su).  A  sample 
stress/strain  curve  for  one  of  the  prism  tests  is  shown  in  Figure  21. 
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Figure  21.  A  sample  stress/strain  curve  for  axial  testing  of  masonry  prisms. 


Numerical  modeling  of  Mehrabi  masonry  prisms 

Two  masonry  prisms,  one  made  of  hollow  blocks  and  the  other  with  solid 
blocks,  were  modeled  using  the  DIANA  program.  Two-dimensional  model¬ 
ing  and  analysis  were  carried  out  using  the  dimensions  and  geometry  de¬ 
scribed  in  the  previous  section.  For  the  prisms  with  hollow  blocks,  equiva¬ 
lent  thickness  of  the  blocks  was  considered  to  be  1.8  in.,  and  the  thickness 
of  the  mortar  joints  to  be  1.25  in.  For  the  prisms  with  solid  blocks,  the 
thickness  of  blocks  was  considered  to  be  3.625  in.,  and  width  of  mortar 
bed  joints  was  3.5  in. 

Masonry  units  were  modeled  using  rotating  smeared  crack  models,  and 
mortar  joints  were  modeled  using  interface  models.  Two-dimensional 
plain-stress,  four-node  elements  with  four  integration  points  were  used  for 
masonry  units.  Two-dimensional  four-node  interface  elements  with  two 
integration  points  were  used  for  mortar  joints.  Material  parameters  used 
for  modeling  of  concrete  in  hollow  and  solid  blocks  are  shown  in  Table  7. 
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Table  7.  Material  parameters  for  concrete  in  masonry  units,  prism  analysis. 


Parameter 

Parameter  Definition 

Value  for  Flollow  Blocks 

Value  for  Solid  Blocks 

E 

Modulus  of  elasticity 

2  E+6  (psi) 

2  E+6  (psi) 

V 

Poisson’s  ratio 

0.16 

0.16 

ft 

Tensile  strength 

240  (psi) 

230  (psi) 

Gt1 

First  mode  fracture  energy 

0.09  (psi-in.) 

0.09  (psi-in.) 

P 

Shear  retention  factor 

NA 

NA 

f'c 

Compressive  strength 

2,400  (psi) 

2,300  (psi) 

Gfc 

Fracture  energy  in 
compression 

22  (psi-in.) 

22  (psi-in.) 

Tension  Curve 

Shape  of  tensile  stress/strain 
curve 

Exponential 

Exponential 

Compression 

Curve 

Shape  of  compressive 
stress/strain  curve 

Parabolic 

Parabolic 

The  parameters  are  defined  also  in  Table  7.  The  values  of  the  parameters 
were  determined  mostly  by  the  material  test  results  and  infilled  frame 
analysis  calibration  process  reported  by  Mehrabi  et  al.  (1994),  and  by  val¬ 
ues  defined  by  verification  studies  reported  in  Chapter  2  (page  21).  For 
simplicity,  interface  normal  and  shear  stiffnesses  for  solid  and  hollow 
blocks  are  assumed  to  be  the  same,  and  minor  differences  introduced  by 
Mehrabi  et  al.  (1994)  in  their  calibration  process  were  ignored.  The  tensile 
strength  of  concrete  for  which  no  test  result  was  reported  is  assumed  to  be 
10%  of  the  compressive  strength.  The  fracture  energy  in  compression,  Gfc , 
was  increased  to  22  psi-in.  to  provide  strain  at  maximum  strength  and  a 
shape  of  descending  branch  (softening)  in  normal  stress/normal  strain 
curve  similar  to  that  obtained  from  the  test  results. 

Material  parameters  and  their  definitions  for  interface  elements  are  shown 
in  Table  8.  The  values  for  the  parameters  are  mostly  determined  according 
to  those  used  by  Mehrabi  et  al.  (1994)  for  bed  joints  in  modeling  of  the  in¬ 
fill  walls,  and  by  values  defined  by  verification  studies  in  the  first  phase  of 
this  study.  Mehrabi  et  al.  used  shear  stiffness  values  for  bed  joints  in 
analysis  of  infill  walls  that  are  considerably  larger  than  those  used  for  their 
verification  studies  for  individual  mortar  joints.  It  was  their  conclusion 
that  lower  values  calibrated  according  to  single  joint  shear  testing  do  not 
correspond  to  practical  values  due  to  the  fact  that  deformation  of  the  test 
machine  had  distorted  the  joint  deformation  measurements.  The  higher 
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stiffness  values  were  the  result  of  their  calibration  efforts  in  modeling  in¬ 
filled  frames. 


Table  8.  Parameters  defining  interface  model,  prism  analysis. 


Parameter  and  definition 

Parameter  value 

Knn  =  normal  numerical  elastic  stiffness  parameter 

74E+4  (psi/in.) 

Kss  =  tangential  numerical  elastic  stiffness  parameter 

90E+4  (psi/in.) 

ft  =  joint  tensile  strength 

40  (psi) 

Gt1  =  first  mode  fracture  energy 

1.61  (psi-in.) 

Co  =  initial  cohesion 

40  (psi) 

(|)i  =  tangent  of  initial  internal  friction  angle 

0.9 

'P  =  tangent  of  dilatancy  angle 

0.005 

(|>r  =  tangent  of  residual  internal  friction  angle 

0.75 

au  =  confining  stress  over  which  the  dilatancy  will  be  zero 

150  (psi) 

5  =  dilatancy  shear  slip  degradation  coefficient 

2.3 

Gt"  =  shear  mode  fracture  energy 

16.1  (psi-in.) 

fc  =  compressive  strength 

2,100  (psi) 

Cs  =  parameter  controlling  the  shear  stress  contribution  to  failure 

1 

Gfc  =  fracture  energy  in  compression,  interface 

55  (psi-in) 

kp  =  norm  of  plastic  strain  associated  with  peak  compressive 
strength 

0.006 

The  compression  cap  value,  f  c,  was  defined  by  the  compressive  strength  of 
masonry  mortar  reported  by  Mehrabi  et  al.  (1994).  The  fracture  energy  in 
compression,  GfC ,  was  increased  to  55  psi-in.  to  provide  strain  at  maxi¬ 
mum  strength  and  a  gradual  descending  branch  (softening)  in  the  normal 
stress/normal  strain  curve  similar  to  that  obtained  from  test  results. 

All  nodes  at  the  bottom  of  the  lower  masonry  unit  were  restrained  in  every 
direction  and  a  uniform  vertical  displacement  was  applied  at  upper  nodes 
in  the  top  masonry  unit.  The  loaded  nodes  were  restrained  from  horizontal 
movement,  assuming  a  perfect  bond  between  loading  cap  and  masonry 
units  in  the  prism  test.  A  sample  finite  element  mesh  and  stress  distribu¬ 
tion  is  shown  in  Figure  22.  The  stress  distribution  reflects  the  formation  of 
shear  cones  in  the  prisms.  Stress  rise  at  the  upper  bed  joint  shows  initia¬ 
tion  of  the  failure  at  that  joint. 
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Figure  22.  A  sample  finite  element  mesh  and  stress  contour  for  masonry  prism  analyses. 


Figure  23  shows  normal  stress/normal  strain  curves  for  hollow  and  solid 
block  prisms.  The  maximum  strength,  strain  at  maximum  strength,  and 
softening  curves  agree  well  with  the  test  results.  The  normal  stress  is  cal¬ 
culated  based  on  an  equivalent  thickness  of  1.8  in.  and  3.625  in.  for  hollow 
and  solid  block  prisms,  respectively.  It  should  be  noted  that  the  failure  of 
the  prisms  in  the  analysis  was  initiated  and  dominated  by  the  failure  of 
mortar  joints.  This  type  of  failure  is  not  precisely  consistent  with  the  labo¬ 
ratory  tests  in  which  failure  is  dominated  by  conic  shear  failure  of  the  ma¬ 
sonry  units.  The  two-dimensional  model  used  for  the  analysis  is  not  capa¬ 
ble  of  generating  the  confinement  in  the  mortar  required  for  transferring 
the  failure  to  the  masonry  blocks.  This  inability  to  confine  the  mortar  is  a 
limitation  of  the  two-dimensional  model  used  in  this  study. 
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Normal  Strain  (in/in) 


Figure  23.  Normal  stress/normal  strain  curves  for  hollow  and  solid  block  masonry  prisms. 

Comparison  of  experimental  and  analytical  results 

A  prism  built  of  three  blocks  and  two  bed  joints  was  tested  in  compression. 
The  blocks  were  modeled  with  the  same  material  parameters  as  above.  The 
joints  were  modeled  using  the  aforementioned  interface  model.  Loading 
was  applied  monotonically.  The  small  relative  stiffness  of  the  joint  with 
respect  to  the  block,  which  could  have  potentially  caused  numerical  prob¬ 
lems,  was  increased  by  a  factor  of  io.  The  material  parameters  for  the 
prism  masonry  assembly,  including  CMU  and  interface  elements  (mortar), 
are  given  in  Table  g.  Stress  contours  are  shown  in  Figure  24  and  Figure  25. 


Table  9.  Material  para  meters  for  the  prism  masonry  assembly. 


CMU  Element  Material  Parameters 

E 

V 

ft 

Gt* 

Tensile 

curve 

Shear* 

curve 

P 

f'c 

Comp. 

Curve 

Gc 

Gt" 

fc 

Cs 

Gfc 

Kp 

3.5E+6 

0.16 

200 

0.03 

Exponential 

BEDIAG 

0.02 

3,000 

Parabola 

22 

16.1 

1,500 

1.0 

2 

0.006 

Interface  Element  Material  Parameters 

Knn 

Kss 

ft 

Gt' 

c 

<t>i 

Yi 

<|>r 

au 

5 

Gt" 

fc 

Cs 

Gtc 

Kp 

2.8E+4 

3.5E+4 

40 

1.69 

26 

0.9 

0.005 

0.8 

150 

0.5 

16.1 

1,500 

1.0 

2 

0.006 

*  Shear  is  not  defined  for  a  rotating  crack  model. 
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iDIANA  8.1-460  :  U.S.Army  Construction  Eng. Res.  2-JUN-2004  02:27  prism3-Step-l 

MODEL:  DIANA-W 

LCl:  LOAD  CASE  1 

step:  1  load:  .ie-2 

ELEMENT  PRINC  STRESS  PMIN 

CALCULATED  FROM:  EL.SXX.G 

MAX  =  -147  MIN  =  -195 

1-151 

1-156 

1-160 

1-164 

1-169 

>-173 

-177 

-182 

-186 

-190 

L 

1 

Figure  24.  The  prism  at  step  1. 


iDIANA  8.1-460  :  U.S.Army  Construction  Eng. Res.  2-JUN-2004  04:16  prism3-peak 

model:  diana-W 

DEFORMATION  =  10 

LCl:  LOAD  CASE  1 

STEP:  23  LOAD:  . 23E-1 

ELEMENT  EL.SXX.G  SYY 

MAX  =  -. 155E4 

MIN  =  -. 207E4 

I-.16E4 

I-.161E4 

■-.161E4 

I-.162E4 

I-.163E4 

I-.163E4 

I-.164E4 

■-.165E4 

-  .166E4 

-  .166E4 

-  .167E4 
.168E4 

-  .168E4 
0-.169E4 

-  .17E4 

-  .17E4 
Q-.171E4 

-  .172E4 

-  .172E4 

-  .173E4 

-  .174E4 
- .174E4 

-  .175E4 
- .176E4 
- .177E4 
- .177E4 

I-.178E4 

I-.179E4 

I-.179E4 

I-.18E4 

L 

1 

Figure  25.  The  prism  at  the  peak  load  (step  23). 


The  diagonal  localization  of  strains  can  be  observed.  These  localization 
bands  unload  faster  than  adjacent  material,  causing  overloading  and  a 
higher  damage  rate  in  the  neighborhood.  Figure  26  shows  successful  cap¬ 
ture  of  peak  and  post  peak  regions. 
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iDIANA  8.1-460  :  U.S.Army  Construction  Eng. Res. 

2-JUN-2O04  04:06  prism3-LD 

model: 

DIANA-W 

MAX/MIN 

on  whole  graph: 

YMAX  = 

787  YMIN  =  53.8 

XMAX  = 

-. IE-2  XMIN  =  -. 5E-1 

VARIATION  OVER  LOADCASES 

NODE  447 
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NODAL  TDTX. 
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Figure  26.  Load-displacement  of  node  447. 

This  behavior  is  equivalent  to  a  peak  of  i,6oo  psi  at  a  strain  of  0.002.  The 
strength  appears  to  be  acceptable,  and  it  is  an  indication  that  the  mortar 
joint  has  affected  the  prism’s  overall  strength.  In  reality,  a  3,000  psi  block 
placed  in  a  prism  assembly  will  exhibit  decreases  in  overall  strength  to 
1,500  psi.  The  peak  load  strain  of  0.002,  however,  is  rather  small,  indicat¬ 
ing  excessive  damage  propagation  within  the  materials.  Experimental  data 
on  different  prisms  indicate  that  this  strain  should  range  between  0.003  ~ 
0.005.  Element  stress  distribution  is  shown  in  Figure  27. 
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iDIANA  8.1-460  :  U.S.Arniy  Construction  Eng. Res.  2-JUN-2004  02:43  prism-3-step6 

model:  diana_w 

LCl :  LOAD  CASE  1 

STEP:  6  LOAD:  . 6E-2 

ELEMENT  PRINC  STRESS  PMIN 

CALCULATED  FROM:  EL.SXX.G 

MAX  =  -789 

MIN  =  -. 105E4 

i  i  i  i  i  i  i  i  i  i 

SSmSmSSSSjiM 

Y 

* 

1 

Figure  27.  Element  stress  distribution  step  6. 
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5  Modeling  and  Analysis  of  Masonry-Infilled 
R/C  Frames 

This  chapter  examines  the  capabilities  and  identifies  the  limitations  of  the 
DIANA  finite  element  program  in  simulating  the  behavior  of  masonry- 
infilled  structures  subjected  to  lateral  loading.  The  purpose  of  the  investi¬ 
gation  is  to  establish  a  framework  for  future  modeling  and  analytical  stud¬ 
ies  of  these  highly  nonlinear  structures.  Two  specimens  among  those 
tested  by  Mehrabi  et  al.  (1994)  with  distinctively  different  load-carrying 
and  failure  mechanisms  are  considered.  Both  specimens  are  based  on 
frames  that  are  not  designed  for  high  seismicity  —  so-called  weak  frames 
—  which  are  susceptible  to  developing  shear  failure  in  the  columns.  This 
type  of  failure  is  of  the  greatest  concern  for  R/C-infilled  frame  structures. 

DIANA  was  first  used  to  model  and  analyze  one  frame  with  a  relatively 
weak  infill  in  order  to  further  calibrate  the  models  for  obtaining  agreement 
between  the  analytical  and  experimental  results.  The  failure  of  such  a 
specimen  is  expected  to  be  governed  by  shearing  and  slip  along  masonry 
bed  joints  (mechanism  5  in  Figure  1).  Then,  the  calibrated  models  were 
used  in  DIANA  to  analyze  the  performance  of  a  second  frame,  this  latter 
specimen  having  a  relatively  strong  infill.  The  failure  mechanism  for  this 
second  frame  would  be  expected  to  be  governed  by  diagonal  cracking  of 
the  infill  and  shear  failure  of  columns  (mechanism  2  in  Figure  1).  The  goal 
of  the  investigation  was  to  verily  that  the  models  calibrated  with  the  re¬ 
sults  of  one  test  could  be  applied  to  a  different  frame  with  repeatable  Valid¬ 
ity- 

Modeling  assumptions 

The  specimens  tested  by  Mehrabi  et  al.  (1994)  and  used  in  this  analytical 
study  were  half-scale  frame  models  representing  the  interior  bay  at  the 
bottom  story  of  a  prototype  frame.  The  prototype  frame  was  a  six-story, 
three-bay,  moment-resisting  R/C  frame  with  a  45  x  15  ft  tributary  floor 
area  at  each  story.  Two  types  of  frames  were  designed  for  the  prototype 
structure  with  respect  to  the  lateral  loadings:  one  with  a  weak  frame  and 
weak  infill,  and  one  with  a  weak  frame  and  strong  infill.  The  weak  infill  de¬ 
sign  was  based  on  the  assumption  of  a  strong  wind  load,  and  the  strong 
infill  design  based  on  the  assumption  of  equivalent  static  forces  for  strong 
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seismic  loading.  The  infill  panels  were  constructed  of  4  x  4  x  8  in.  (nomi¬ 
nal)  hollow-core  and  solid  concrete  masonry  blocks  to  represent  the  weak 
and  strong  infill  panels,  respectively. 

The  accuracy  and  reliability  of  an  analytical  model  for  simulating  the  be¬ 
havior  of  an  infilled  frame  strongly  depend  on  the  capability  of  the  model 
to  predict  the  load-carrying  and  failure  mechanisms  as  well  as  estimate  the 
strength  and  deformations.  The  failure  mechanisms  of  these  specimens 
correspond  with  mechanisms  5  and  2  in  Figure  1,  respectively.  Those 
mechanisms  are  common  failures  for  frames  with  unreinforced  masonry 
infill  that  are  not  designed  for  high  seismicity,  in  accordance  with  recent 
design  codes. 

Geometry  and  details  of  the  selected  specimens  are  shown  in  Figure  28. 
Each  test  specimen  was  subjected  to  constant  vertical  loading  and  mono- 
tonically  increasing  lateral  loading  with  the  scheme  shown  in  Figure  29.  In 
this  figure,  P2  is  22  kips,  P3  is  11  kips,  and  d  is  equal  to  16.5  in.  Figure  30 
and  Figure  31  show  damage  to  the  two  specimens  mapped  after  comple¬ 
tion  of  the  tests  for  weak  and  strong  infill  specimens,  corresponding  to 
mechanisms  5  and  2,  respectively,  in  Figure  1. 


Figure  28.  Geometry  and  details  of  test  specimen  (Mehrabi  et  al.  1994). 
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Figure  30.  Failure  pattern  from  laboratory  test  of  frame  with  weak  infill  (Mehrabi  et  al.  1994). 
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Figure  31.  Failure  pattern  from  laboratory  test  of  frame  with  strong  infill  (Mehrabi  et  al.  1994). 

Figure  32  and  Figure  33  show  deformed  shapes  of  the  finite  element  mod¬ 
els  for  the  weak  infill  specimen  and  strong  infill  specimen,  respectively. 
These  figures  show  opening  and  slip  along  the  mortar  joints  for  the  weak 
infill  specimen,  and  diagonal  cracking  and  opening  of  mortar  joints  in  in¬ 
fill  and  shear  failure  of  the  top  end  of  the  windward  column  for  the  speci¬ 
men  with  strong  infill.  The  numerical  responses  for  the  simulations  were 
in  good  agreement  with  experimental  test  results. 
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Figure  32.  Deformed  shape  predicted  by  analysis  of  frame 
with  weak  infill  (Mehrabi  et  al.  1994). 


Figure  33.  Deformed  shape  predicted  by  analysis  of  frame 
with  strong  infill  (Mehrabi  et  al.  1994). 
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Development  of  infill  panel  model 

Frame  and  infill  elements 

The  infilled  frame  with  geometry  and  details  shown  in  Figure  28  was  mod¬ 
eled  with  the  DIANA  finite  element  program.  Two-dimensional,  plain- 
stress,  four-node  elements  with  four  integration  points  were  used  for 
modeling  the  concrete  in  R/C  frame  and  masonry  units.  Two-dimensional, 
four-node  interface  elements  with  two  integration  points  were  used  for 
modeling  mortar  bed  and  head  joints  and  joints  between  the  infill  and  the 
frame.  DIANA  rotating-smeared-crack  constitutive  models  for  modeling 
the  concrete  and  DIANA  interface  elements  for  modeling  the  mortar  joints 
were  utilized,  as  described  previously.  Reinforcement  bars  for  the  frame 
were  modeled  with  DIANA  elastic-hardening  plastic,  two-node  discrete 
bar  elements.  The  loading  plates  were  modeled  with  linear  elastic,  four- 
node,  plane-stress  elements.  The  top  loading  beam  was  modeled  with  two- 
node  beam  elements  using  linear-elastic  material  with  steel  properties, 
while  the  frame  footing  was  modeled  with  four-node  plane  stress  elements 
with  linear-elastic  material  with  concrete  properties.  The  finite  element 
mesh  is  shown  in  Figure  34  at  initial  stages  of  lateral  loading.  Each  ma¬ 
sonry  unit  was  discretized  into  two  elements  with  an  aspect  ratio  of  1.  The 
loading  scheme  followed  that  illustrated  in  Figure  29.  Vertical  load  was 
first  applied  in  increments  to  the  maximum  and  then  kept  constant  while 
lateral  displacement  loading  was  applied  gradually. 

Parameter  setting  and  analysis 

When  analyzing  infilled  frames,  Mehrabi  et  al.  (1994)  adjusted  the  normal 
and  shear  stiffnesses  of  mortar  joints  to  about  30  times  those  calibrated  by 
their  laboratory  direct  shear  test  results.  The  need  for  adjustment  was  at¬ 
tributed  to  the  inaccuracy  of  interface  normal  and  shear  displacements  in 
the  elastic  region  of  laboratory  test  responses  that  occurred  because  of  de¬ 
flection  of  the  test  fixture.  They  differentiated  among  mortar  bed  joints, 
head  joints,  and  joints  between  the  wall  and  frame  by  introducing  different 
interface  material  parameters  and  thicknesses. 
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Parameter  setting  for  frame  with  weak  infill 

For  the  first  trial  using  DIANA,  material  parameters  were  selected  to  agree 
with  those  used  by  Mehrabi  et  al.  (1994).  The  parameters  for  the  bed  joints 
were  those  used  for  the  prism  analysis  previously  described,  also  reflected 
in  Table  3.  The  use  of  the  normal  and  shear  stiffnesses  from  Table  10  re¬ 
sulted  in  divergence  and  a  lack  of  solution  at  initial  stages  of  the  analysis. 
This  divergence  occurred  with  the  initiation  of  shear  cracking  and  slip  at 
interfaces.  The  numerical  process  of  determining  the  stress  on  the  yield 
surface  (return  mapping)  failed  in  tension-shear  or  compression-shear 
corner  zones.  Persistent  efforts  to  prevent  the  divergence  of  the  algorithm 
with  reduction  of  step  sizes  and  tolerances  in  the  practical  range,  and  ap¬ 
plication  of  various  available  solution  methods  did  not  resolve  the  conver¬ 
gence  issue.  It  was  concluded  that  sharp  corners,  especially  at  tension  cut¬ 
off  zones,  need  to  be  repaired  in  the  yield  surface  of  the  interface  model  if 
this  problem  is  to  be  avoided.  Mehrabi  et  al.  (1994)  utilized  a  hyperbolic 
yield  surface  which  avoided  corners  in  the  shear-tension  zone,  and  their 
model  did  not  include  a  compression  cap. 
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Table  10.  Material  parameters  for  concrete  in  frame  and  masonry  units. 


Parameter 

Description 

Concrete  in 

Frame 

Concrete  in 
masonry  hollow 
units 

Concrete  in 
masonry  solid 
units 

E 

Modulus  of  elasticity 

3.55  E+6  (psi) 

2.0  E+6  (psi) 

2.0  E+6  (psi) 

V 

Poisson’s  ratio 

0.16 

0.16 

0.16 

ft 

Tensile  strength 

390  (psi) 

240  (psi) 

230  (psi) 

Gf' 

First  mode  fracture  energy 

0.09  (psi-in) 

0.09  (psi-in) 

0.09  (psi-in) 

P 

Shear  retention  factor 

NA 

NA 

NA 

f'c 

Compressive  strength 

3,900  (psi) 

2,400  (psi) 

2,300  (psi) 

Gtc 

Fracture  energy  in 
compression 

22  (psi-in) 

22  (psi-in) 

22  (psi-in) 

Tension  Curve 

Shape  of  tensile 
stress/strain  curve 

Exponential 

Exponential 

Exponential 

Compression  Curve 

Shape  of  compressive 
stress/strain  curve 

Parabolic 

Parabolic 

Parabolic 

However,  during  verification  of  constitutive  models  in  the  first  phase  of 
this  study,  with  the  set  of  parameters  calibrated  on  the  basis  of  the  labora¬ 
tory  material  test  results,  the  analysis  did  not  encounter  a  convergence 
problem.  That  result  was  due  to  much  lower  normal  and  shear  stiffness 
values  than  those  utilized  by  Mehrabi  et  al.  (1994)  for  their  infilled  frame 
analyses.  Therefore,  it  was  decided  to  carry  out  the  analysis  with  lower 
stiffness  values  according  to  those  used  in  the  material-level  investigation. 
Another  trial  analysis  using  DIANA  confirmed  conclusions  by  Mehrabi  et 
al.  (1994)  that  the  use  of  these  low  stiffness  values  does  not  yield  an 
agreement  between  experimental  and  analytical  responses.  This  trial 
analysis  clearly  shows  that  low  normal  stiffness  for  interfaces  resulted  in 
transferring  much  of  the  vertical  load  to  the  frame  columns  and  less  of  it  to 
the  infill  wall.  Since  the  shear  resistance  of  the  infill  wall  strongly  depends 
on  normal  stresses,  low  normal  stresses  on  infill  resulted  in  significantly 
lower  lateral  resistance  for  the  infilled  frame.  Also,  lower  shear  stiffness 
for  infill  mortar  joints  resulted  in  a  significantly  lower  initial  stiffness  for 
the  analytical  response  curve  than  seen  in  the  experimental  results.  To 
find  an  agreement  between  analytical  and  experimental  results,  higher 
stiffnesses  need  to  be  used. 

To  use  a  higher  stiffness  for  interfaces  and  to  avoid  the  numerical  conver¬ 
gence  problem,  researchers  decided  to  keep  the  stiffnesses  at  the  lower 
values  for  all  interfaces  as  the  base  stiffness  parameters  (for  which  a  good 
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numerical  convergence  is  guaranteed  at  material  level)  and  to  increase  the 
theoretical  width  of  joints  to  provide  the  higher  overall  stiffness  required 
for  obtaining  agreement  with  the  experimental  results.  However,  to  avoid 
an  unwanted  and  unrealistic  increase  in  the  overall  interface  strength  val¬ 
ues  due  to  increase  in  width,  the  parameters  affecting  the  strength  (i.e., 
compressive  strength,  tensile  strength,  and  cohesion)  have  to  be  reduced 
by  the  same  proportion  the  thicknesses  are  increased.  The  confining 
stress,  ctu,  is  related  only  to  evolution  of  dilatancy,  so  it  is  not  reduced. 

After  a  series  of  trial  runs  with  various  multiples  of  width  increase,  a  ten¬ 
fold  increase  in  the  stiffness  of  bed  joints  over  the  base  values  (estimated 
originally  by  calibration  using  material  test  results)  was  found  to  provide  a 
good  agreement  between  the  analytical  and  experimental  responses  for  the 
frame  with  weak  infill.  The  stiffness  of  the  head  joints  and  the  joints  be¬ 
tween  the  frame  and  the  infill  were  adjusted  accordingly,  keeping  the  same 
proportion  considered  by  Mehrabi  et  al.  (1994)  in  their  modeling  effort 
with  respect  to  the  bed  joints.  Material  parameters  for  the  concrete  mate¬ 
rial  in  the  frame  and  masonry  units  were  calculated  according  to  the  test 
results  provided  by  Mehrabi  et  al.  (1994)  and  the  verification  process  in 
this  study.  Table  10  shows  material  parameters  and  their  definition  for  the 
concrete  used  in  the  frame  and  the  masonry  units.  Table  11  shows  the  tar¬ 
get  material  parameters  and  the  parameters  actually  used  for  bed  joints, 
head  joints,  and  joints  between  frame  and  wall,  respectively. 


Table  11.  Material  parameters  for  hollow  block  infill,  DIANA  interface  model. 


Parameter  Set 

Width 

(in.) 

Knn 

(psi) 

Kss 

(psi) 

ft 

(psi) 

Co 

(psi) 

tg<t* 

tgf 

f'c 

(psi) 

Target  values  for  bed 
joints 

1.25 

280,000 

350,000 

40 

40 

0.9 

0.75 

1,500 

Actual  values  used 
for  bed  joints 

12.5 

28,000 

35,000 

4 

4 

0.9 

0.75 

150 

Target  values  for 
head  joint 

1.25 

215,300 

269,200 

10 

10 

0.8 

0.7 

1,500 

Actual  values  used 
for  head  joint 

12.5 

21,530 

26,920 

1 

1 

0.8 

0.7 

150 

Target  values 

For  frame/ 
wall  joints 

1.4 

215,300 

269,200 

20 

20 

0.8 

0.7 

1,500 

Actual  values  used 
for  frame/wall  joints 

14 

21,530 

26,920 

2 

2 

0.8 

0.7 

150 
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Analysis  results  for  frame  with  weak  infill 

Within  the  geometry  and  material  properties  described  above,  the  analysis 
of  the  selected  masonry-infilled  R/C  frame  with  weak  infill  was  carried  out 
for  a  maximum  lateral  displacement  of  0.8  in.  Figure  34  shows  the  princi¬ 
pal  stress  contour  at  0.07  in.  lateral  displacement,  indicating  signs  of 
forming  the  diagonal  stmts  that  are  the  load  carrying  mechanism  of  in¬ 
filled  frames  at  initial  loading  stages. 


Figure  35  shows  the  deformed  shape  and  principal  stress  contour  at  0.31 
in.  lateral  deflection  before  maximum  strength  is  reached.  This  figure 
shows  separation  and  slip  at  bed  and  head  joints,  which  is  the  governing 
failure  mechanism  for  this  infilled  frame.  It  also  shows  that  stress  at 
loaded  corners  of  the  infill  has  reached  nearly  the  compressive  strength  of 
the  masonry,  signaling  crushing  at  loaded  corners  under  higher  lateral 
displacements.  These  results  agree  well  with  the  experimental  and  analyti¬ 
cal  results  by  Mehrabi  et  al.  (1994)  and  show  the  capability  of  the  models 
in  predicting  the  load  carrying  and  failure  mechanisms  of  the  R/C  frame 
with  weak  masonry  infill. 
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Figure  35.  Deformed  shape  and  stress  contour 
for  frame  with  weak  infill,  analysis  using  DIANA. 


Figure  36  shows  response  curves  obtained  from  this  analysis  compared 
with  those  from  the  laboratory  test.  Stiffness,  response  trend  (including 
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initiation  of  nonlinear  behavior  due  to  mortar  joint  shear  cracking  and 
separation),  and  lateral  resistance  obtained  numerically  agree  well  with 
the  experimental  results. 


Figure  36.  Experimental  and  numerical  (DIANA)  lateral  load/lateral  displacement  curves 

for  frame  with  weak  infill. 

Parameter  setting  for  frame  with  strong  infill 

To  examine  the  capability  of  the  calibrated  models  to  predict  the  behavior 
of  an  infilled  frame  with  different  characteristics,  the  FE  model  created 
and  calibrated  for  analysis  of  the  frame  with  weak  infill  was  used  for  analy¬ 
sis  of  the  frame  with  strong  infill.  The  only  parameter  in  the  interface 
model  needing  adjustment  was  the  compressive  strength  (or  the  compres¬ 
sion  cap),  f  c,  to  reflect  the  higher  strength  of  a  masonry  assembly  made  of 
solid  blocks.  The  width  of  mortar  joints  also  had  to  be  adjusted  to  higher 
values.  These  changes  are  shown  in  Table  12.  A  thickness  of  3.625  in.  was 
used  for  masonry  units  in  the  infill.  The  model  was  then  analyzed  using 
the  same  loading  scheme  that  was  applied  to  the  frame  with  weak  infill. 


Table  12.  Material  parameters  for  mortar  joints  in  solid  block  infill 
that  are  different  from  hollow  block  infill,  DIANA  interface. 


Parameter 

Set 

Width  of  Bed 
Joints  (in.) 

Width  of  Head 

Joints  (in.) 

Width  of  Frame  to 
Wall  Joints  (in.) 

f'c  (psi) 

Target  values 

3.5 

3 

3.5 

2,000 

Actual  values  used 

35 

30 

35 

200 
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Analysis  results  for  frame  with  strong  infill 

The  analysis  results  indicated  that  the  predicted  behavior  of  the  frame 
with  strong  infill  using  the  models  described  above  agreed  well  with  the 
experimental  behavior  for  the  initial  portion  of  the  response  curve,  i.e.,  for 
lateral  displacement  smaller  than  0.2  in.  The  analytical  response  did  not 
flatten,  however,  and  the  lateral  resistance  continued  to  increase  with  in¬ 
creasing  displacement  far  beyond  the  strength  obtained  in  the  experimen¬ 
tal  test  (Figure  37).  Review  of  the  numerical  results  indicated  that  the 
shear  failure  at  the  top  end  of  the  windward  column,  expected  to  govern 
the  failure  of  the  infilled  frame  based  on  the  experimental  results,  never 
occurred.  Furthermore,  with  stronger  mortar  joints  in  strong  infill,  sliding 
of  bed  joints  did  not  occur.  Only  some  diagonal  cracking  and  mortar  joint 
separation  were  observed.  This  behavior  resulted  in  development  of  a  con¬ 
fined  diagonal  compression  stmt  in  the  infill,  which  provided  much  higher 
lateral  strength  than  expected.  This  problem  was  attributed  to  the  ineffi¬ 
ciency  of  the  smeared  crack  model  for  concrete  in  the  R/C  frame  for  mod¬ 
eling  the  shear  failure. 


Figure  37.  Experimental  and  numerical  (DIANA)  lateral  load/lateral  displacement  curves 

for  frame  with  strong  infill. 
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To  overcome  this  modeling  problem,  a  discontinuity  in  the  form  of  an  in¬ 
terface  model  was  introduced  at  the  column  ends.  This  allowed  the  shear 
failure  of  the  columns  along  these  interfaces,  as  it  was  anticipated  for  this 
infilled  frame.  The  material  properties  of  the  interface  element  at  the  col¬ 
umn  ends  were  similar  to  those  of  the  interfaces  in  the  masonry.  However, 
some  parameters,  e.g.,  tensile  strength,  cohesion,  and  compressive 
strength,  were  adjusted  to  reflect  properties  of  concrete  in  the  frame.  To 
avoid  convergence  problems,  the  method  of  increasing  interface  thick¬ 
nesses  adopted  for  the  case  of  interfaces  in  infill  was  used  here  as  well. 


With  the  above  modification,  the  model  was  analyzed  again.  As  it  can  be 
seen  in  Figure  37,  this  time  the  response  agreed  well  with  the  experimental 
results.  Figure  38  shows  the  deformed  shape  and  principal  stress  contour 
at  0.215  in.  lateral  displacement  before  shear  failure  of  the  column  occurs. 
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LC2:  LOAD  CASE  2 
STEP:  265  LOAD:  .215 
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Figure  38.  Deformed  shape  and  stress  contour  for  frame 
with  strong  infill  at  0.215  in.  lateral  displacement,  analysis  using  DIANA. 


Separation  of  head  and  bed  joints  in  a  diagonal  pattern  and  development 
of  compression  zones  at  loaded  corners  can  be  seen  in  Figure  38.  Figure  39 
shows  the  deformed  shape  and  displacement  field  contour  at  0.795  in-  lat¬ 
eral  displacement.  This  figure  shows  the  shear  failure  of  the  top  end  of  the 
windward  column,  as  illustrated  also  by  the  experimental  results  in  Figure 
31.  A  large  part  of  the  lateral  displacement  has  been  absorbed  by  the  shear 
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failure  of  the  column  and  displacement  of  the  upper  right  portion  of  the 
frame  and  the  infill.  Load  carrying  and  failure  mechanisms  of  the  modified 
model  agree  well  with  the  experimental  results. 


MODEL:  WEAKPARAM_W 
DEFORMATION  =  5 
LC2:  LOAD  CASE  2 
STEP:  845  LOAD:  .795 
NODAL  TDTX. . .G  RESTDT 
MAX  =  .93  MIN  =  0 
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Figure  39.  Deformed  shape  and  displacement  field  contour  for  frame  with  strong  infill  at 
0.795  in.  lateral  displacement,  analysis  using  DIANA. 
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6  Conclusions 

The  objective  of  this  work  was  to  investigate  the  utility  and  effectiveness  of 
numerical  finite  element  modeling,  simulation,  and  analysis  to  predict  the 
performance  characteristics  of  masonry-infilled  wall  panels  and  constitu¬ 
ent  materials  subjected  to  lateral  loading.  One  purpose  of  the  work  was  to 
determine  whether  numerical  analysis,  based  on  well  calibrated  finite  ele¬ 
ment  models  and  parametric  studies,  can  serve  as  a  valid  substitute  for  ex¬ 
pensive  physical  modeling  and  testing  of  concrete  and  masonry  structures. 

This  investigation  identified  a  modeling  approach  that  addresses  the  dis¬ 
continuities  in  masonry  by  introducing  interface  elements  that  help  to  rep¬ 
resent  the  masonry  joints  more  accurately  than  in  the  past.  Through  an 
extensive  literature  review,  the  authors  identified  constitutive  material 
models  for  numerical  analysis  of  R/C  frames  infilled  with  unreinforced 
concrete  masonry  blocks,  and  demonstrated  their  capabilities  through 
verification  studies.  A  combination  of  continuum  and  interface  constitu¬ 
tive  models  was  selected  as  the  most  appropriate  approach. 

In  order  to  facilitate  the  use  of  these  models  by  researchers  and  designers, 
a  commercial  finite  element  program  with  compatible  capabilities  was 
identified.  That  program,  DIANA,  was  employed  to  test  the  constitutive 
models  using  data  from  previous  experimental  tests,  and  to  assess  the 
strengths,  weaknesses,  capabilities,  and  limitations  of  those  models. 

First  DIANA  was  used  to  analyze  masonry  prisms  constructed  of  hollow 
and  solid  blocks.  Next,  a  frame  with  weak  infill  (hollow  blocks)  was  mod¬ 
eled  and  analyzed,  and  the  results  were  used  to  perform  further  calibration 
of  the  models  to  obtain  agreement  between  analytical  and  experimental 
results.  The  failure  of  this  modeled  frame  was  governed  by  shear  cracking 
and  slip  along  masonry  bed  joints.  Then  a  second  frame,  one  with  a  strong 
infill  (solid  blocks),  was  developed  using  the  calibrated  models,  and  ana¬ 
lyzed  using  DIANA.  The  goal  was  to  verily  that  the  constitutive  models 
calibrated  with  the  results  of  one  test  would  be  valid  when  used  in  a  differ¬ 
ent  structural  model.  The  failure  mechanism  for  this  second  frame  was  ex¬ 
pected  to  be  governed  by  diagonal  cracking  in  the  infill  and  shear  failure  at 
the  top  end  of  the  windward  column.  However,  modeling  of  the  second 
frame  required  the  introduction  of  interface  elements  in  the  columns  to 
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allow  shear  failure  of  the  columns,  which  did  not  occur  as  expected  using 
only  the  smeared  crack  model  for  the  concrete  material  in  the  frame. 

The  models  used  in  DIANA  showed  good  capabilities  for  modeling  and 
predicting  load-carrying  and  failure  mechanisms  for  such  complex  struc¬ 
tures  as  infilled  frames.  Such  capabilities  are  recognized  to  be  critical  for 
reliably  predicting  the  strength  and  ductility  of  R/C  infilled  frame  designs. 
A  significant  limitation  of  the  models  in  DIANA  is  a  convergence  problem 
with  the  numerical  scheme  in  return  mapping  on  the  yield  surface  for 
higher  stiffnesses  for  interface  elements,  which  is  a  significant  impediment 
for  a  robust  and  reliable  solution.  A  short-term  solution  was  implemented 
to  overcome  this  problem.  It  is  believed  that  the  numerical  inefficiency  of 
the  models  can  be  overcome  by  improvements  to  the  failure  surface  and 
return  mapping  schemes  of  the  constitutive  model  for  interface  elements. 
For  the  case  of  potential  shear  failure  in  columns,  it  was  shown  that  the 
use  of  interface  elements  in  the  column  ends  can  compensate  for  inability 
of  the  smeared  crack  formulation  to  model  shear  failure.  In  any  case,  it  is 
realized  that  for  successful  use  of  any  FE  program,  the  model  parameters 
need  to  be  calibrated  using  appropriate  material-  and  structural-level  ex¬ 
perimental  results.  Once  calibrated,  the  models  can  be  used  reliably  for 
parametric  studies  of  the  behavior  of  infilled  frames  subjected  to  lateral 
loading. 
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simple,  accurate  infill  masonry  model.  The  objective  of  this  study  was  to  identify  suitable  numerical  constitutive  models  and  demon¬ 
strate  their  capabilities.  Two  were  studied  here:  a  cohesive  interface  model  to  simulate  the  behavior  of  mortar  joints  between  masonry 
units  as  well  as  the  behavior  of  the  frame-to-panel  interface;  and  a  smeared  crack  finite  element  formulation.  The  interface  model  was 
able  to  account  for  the  shearing,  residual  shear  strength,  and  opening  and  closing  of  joints  under  cyclic  shear  loads  in  simple  combi¬ 
nations  of  concrete  blocks  and  mortar  joints. 
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